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The dynamical behavior of open quantum systems plays a key role in many applica¬ 
tions of quantum mechanics, examples ranging from fundamental problems, such as the 
environment-induced decay of quantum coherence and relaxation in many-body sys¬ 
tems, to applications in condensed matter theory, quantum transport, quantum chem¬ 
istry and quantum information. In close analogy to a classical Markovian stochastic 
process, the interaction of an open quantum system with a noisy environment is often 
modeled phenomenologically by means of a dynamical semigroup with a corresponding 
time-independent generator in Lindblad form, which describes a memoryless dynamics 
of the open system typically leading to an irreversible loss of characteristic quantum fea¬ 
tures. However, in many applications open systems exhibit pronounced memory effects 
and a revival of genuine quantum properties such as quantum coherence, correlations 
and entanglement. Here, recent theoretical results on the rich non-Markovian quantum 
dynamics of open systems are discussed, paying particular attention to the rigorous 
mathematical definition, to the physical interpretation and classification, as well as to 
the quantification of quantum memory effects. The general theory is illustrated by 
a series of physical examples. The analysis reveals that memory effects of the open 
system dynamics reflect characteristic features of the environment which opens a new 
perspective for applications, namely to exploit a small open system as a quantum probe 
signifying nontrivial features of the environment it is interacting with. This article fur¬ 
ther explores the various physical sources of non-Markovian quantum dynamics, such 
as structured environmental spectral densities, nonlocal correlations between environ¬ 
mental degrees of freedom and correlations in the initial system-environment state, in 
addition to developing schemes for their local detection. Recent experiments addressing 
the detection, quantification and control of non-Markovian quantum dynamics are also 
briefly discussed. 


PACS numbers: 03.65.Yz, 03.65.Ta, 03.67.-a, 42.50.-p, 42.50.Lc 
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I. INTRODUCTION 

The observation and experimental control of charac¬ 
teristic quantum properties of physical systems is of¬ 
ten strongly hindered by the coupling of the system 
to a noisy environment. The unavoidable interaction 
of the quantum system with its surroundings generates 
system-environment correlations leading to an irretriev¬ 
able loss of quantum coherence. Realistic quantum me¬ 
chanical systems are thus open systems governed by a 
non-unitary time development which describes all fea¬ 
tures of irreversible dynamics such as the dissipation 
of energy, the relaxation to a thermal equilibrium or a 


stationary nonequilibrium state and the decay of quan- 


turn coherences and correlations ( 

Alicki and Lendi 

1987 

Breuer and Petruccione 

2002 

Davies 1976 

I- 


There is a well-established treatment of the dynam¬ 
ics of open quantum systems in which the open system’s 
time evolution is represented by a dynamical semigroup 
and a corresponding master equation in Lindblad form 


(Gorini et al ., 1976 Lindblad, 1976). If one adopts a mi¬ 


croscopic system-environment approach to the dynamics 
of open systems, such a master equation may be derived, 
for example, with the help of the Born-Markov approx¬ 
imation by assuming a weak coupling between the sys¬ 
tem and its environment. However, it turns out that for 
many processes in open quantum systems the approxi¬ 
mations underlying this approach are not satished and 
that a description of the dynamics by means of a dy¬ 
namical semigroup fails. Typically, this is due to the fact 
that the relevant environmental correlation times are not 
small compared to the system’s relaxation or decoherence 
time, rendering impossible the standard Markov approx¬ 
imation. The violation of this separation of time scales 
can occur, for example, in the cases of strong system- 
environment couplings, structured or finite reservoirs, 
low temperatures, or large initial system-environment 
correlations. 

If the dynamics of an open quantum system substan¬ 
tially deviates from that of a dynamical semigroup one 
often speaks of a non-Markovian process. This term 
refers to a well-known concept of the theory of classical 
stochastic processes and is used to loosely indicate the 
presence of memory effects in the time evolution of the 
open system. However, the classical notions of Marko- 
vianity and non-Markovianity cannot be transferred to 
the quantum regime in a natural way since they are based 
on a Kolmogorov hierarchy of joint probability distribu- 


tions which does not exist in quantum theory ( 

Accardi 

et al. 

1982 

Lindblad 1979 

Vacchini et al. 

20111 

There- 


fore, the concept of a quantum non-Markovian process 
requires a precise definition which cannot be based on 
classical notions only. Many important questions need 
to be discussed in this context: How can one rigorously 
define non-Markovian dynamics in the quantum case, 
how do quantum memory effects manifest themselves in 
the behavior of open systems, and can such effects be 
uniquely identified experimentally through monitoring of 
the open system? The definition of non-Markovianity 
should provide a general mathematical characterization 
which does not rely on any specific representation or ap¬ 
proximation of the dynamics, e.g. in terms of a pertur¬ 
bative master equation. Moreover, definitions of non- 
Markovianity should lead to quantitative measures for 
the degree of non-Markovianity, enabling the compari¬ 
son of the amount of memory effects in different physical 
systems. 

Recently, a series of different proposals to answer these 
questions has been published, rigorously defining the bor- 
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der between the regions of Markovian and non-Markovian 
quantum dynamics and developing quantitative measures 
for the degree of memory effects (see, e.g., (jBreuer et al. 
2009[ Chrusciiiski et alj 2011| [Chruscihski and Manis- 


calco[ 


2014 Luo et al. , 2012[ Rivas et al.[ 2010[ Wolf et al. 


2008), the tutorial paper (Breuer, 2012) and the recent 


review by Rivas et al. (2014)). Here, we describe and 
discuss several of these ideas, paying particular attention 
to those concepts which are based on the exchange of in¬ 
formation between the open system and its environment, 
and on the divisibility of the dynamical map describing 
the open system’s time evolution. We will also explain 
the relations between the classical and the quantum no¬ 
tions of non-Markovianity and develop a general classi¬ 
fication of quantum dynamical maps which is based on 
these concepts. 

The general theory will be illustrated by a series of 
examples. We start with simple prototypical models 
describing pure decoherence dynamics, dissipative pro¬ 
cesses, relaxation through multiple decay channels and 
the spin-boson problem. We then continue with the 
study of the dynamics of open systems which are coupled 
to interacting many-body environments. The examples 
include an Ising and a Heisenberg chain in transverse 
magnetic fields, as well as an impurity atom in a Bose- 


Einstein condensate (Apollaro et al. 2011 Haikka et al. 


2012, 2011). The discussion will demonstrate, in partic¬ 


ular, that memory effects of the open system dynamics 
reflect characteristic properties of the environment. This 
fact opens a new perspective, namely to exploit a small 
open system as a quantum probe signifying nontrivial fea¬ 
tures of a complex environment, for example the critical 


point of a phase transition (Gessner et al. , 2014a Smirne 


et al. 2013). Another example to be discussed here is 


the use of non-Markovian dynamics to determine nonlo¬ 
cal correlations within a composite environment, carrying 
out only measurements on the open system functioning 


as quantum probe (Laine et al. 

2012 

Liu et al. 

2013a 

WiBmann and Breuer 

2014 

)• 


A large variety of further applications of quantum 
memory effects is described in the literature. Interested 
readers may find examples dealing with, e.g., phenomeno¬ 


logical master equations (Mazzola et al. 2010), chaotic 


systems (Znidaric et al.. 2011), energy transfer processes 


in photosynthetic complexes (Rebentrost and Aspuru- 


Guzik 2011), continous variable quantum key distribu- 
tion ( |Vasile et al. 2011b), metrology (Ghin et a/. 2012), 
steady state entanglement (Huelga et al. , 2012), Goulomb 


(Haikka et al. 2013). 


crystals (Borrelli et at., 2013), symmetry breaking (Ghan- 


cellor et al. 2013), and time-invariant quantum discord 


The standard description of the open system dynamics 
in terms of a dynamical map is based on the assumption 
of an initially factorizing system-environment state. The 
approach developed here also allows to investigate the 
impact of correlations in the initial system-environment 


state and leads to schemes for the local detection of 
such correlations through monitoring of the open system 


Gessner and Breuer 

2011 

Laine et al. 

2010b 


In recent years, several of the above features of non- 
Markovian quantum dynamics have been observed exper¬ 


imentally in photonic and trapped ion systems (Gialdi 


et al.\ 

M4) 

Gessner et al.\ 2014b 

Li et al. 

2011 

Liu 

et al. 

2011 

Smirne et al. 

2011 

Tang et al. 

2014D. 

We 


briefly discuss the results of these experiments which 
demonstrate the transition from Markovian to non- 
Markovian quantum dynamics, the non-Markovian be¬ 
havior induced by nonlocal environmental correlations, 
and the local scheme for the detection of nonclassical ini¬ 
tial system-environment correlations. 


II. DEFINITIONS AND MEASURES FOR QUANTUM 
NON-MARKOVIANITY 

A. Basic concepts 

1. Open quantum systems and dynamical maps 


An open quantum system S 

(Alicki and Lendi 

1987 

Breuer and Petruccione 

2002 

Davies 

1976) can 

3e re- 


garded as subsystem of some larger system composed of 
S and another subsystem E, its environment, see Fig. [l] 
The Hilbert space of the total system S + E is given by 
the tensor product space 

T-Lse = T-Ls®T-Le, ( 1 ) 

where Tis and 'He denote the Hilbert spaces of S and 
A, respectively. Physical states of the total system are 
represented by positive trace class operators psE on 'Hse 
with unit trace, satisfying psE ^ 0 and tT:psE = 1- Given 
a state of the total system, the corresponding states of 
subsystems S and E are obtained by partial traces over 
'He and 'Hs, respectively, i.e., we have ps = t^EPSE and 
Pe = trspsE). We denote the convex set of physical states 
belonging to some Hilbert space 'H hy S{'H). 





environment 

interaction Hj 

Hilbert space 'H.e ^ 

open system 


Hilbert space Hs 


Hamiltonian He ^ 

Hamiltonian Hs 


environment state pE 

system state ps 


1. 


> 


FIG. 1 (Color online) Sketch of an open quantum system 
described by the Hilbert space 'Hs and the Hamiltonian Hs, 
which is coupled to an environment with Hilbert space He 
and Hamiltonian He through an interaction Hamiltonian Hi. 

We suppose that the total system S + E is closed and 
governed by a Hamiltonian of the general form 

H = Hs 'S' Ie + Is He + Hj, 


( 2 ) 
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where Hs and He are the free Hamiltonians of system 
and environment, respectively, and H: is an interaction 
Hamiltonian. The corresponding unitary time evolution 
operator is thus given by 

U (t) = exp{—iHt) (^h=l). (3) 

The dynamics of the total system states is obtained from 
the von Neumann equation, 

~j^PSE{t) = -i[H,psE{t)], (4) 

which yields the formal solution 

PSE{t) =UtPSE{0) = U{t)psE{0)U^t). (5) 

An important concept in the theory of open quantum 
systems is that of a dynamical map. To explain this con¬ 
cept we assume that the initial state of the total system 
is an uncorrelated tensor product state 

PS£(0) = Ps(0) ® Pe(0), (6) 


which leads to the following expression for the reduced 
open system state at any time t ^ 0, 

ps{t) = He {C/(t)ps(0) (g)pij(0)C/'l’(t)} . (7) 

Considering a fixed initial environmental state p_e( 0 ) and 
any fixed time t ^ 0, Eq. 0 defines a linear map 

: SiHs) S{Hs) (8) 

on the open system’s state space SiHs) which maps any 
initial open system state ps(0) to the corresponding open 
system state ps(t) at time t: 

Ps{0) Psit) = ^tps{0)- ( 9 ) 


is called quantum dynamical map. It is easy to verify 
that it preserves the Hermiticity and the trace of oper¬ 
ators, and that it is a positive map, i.e., that it maps 
positive operators to positive operators. Thus, maps 
physical states to physical states. 

A further important property of the dynamical map 
is that it is not only positive but also completely pos¬ 
itive. Maps with this property are also known as trace 
preserving quantum operations or quantum channels in 
quantum information theory. Let us recall that a linear 
map is completely positive if and only if it admits a 
Kraus representation (Kraus 1983), which means that 
there are operators Hi on the underlying Hilbert space 
Tis such that $A = that the condition 

of trace preservation takes the form = Is- An 

equivalent definition of complete positivity is the follow¬ 
ing. We consider for any number n = 1,2,... the tensor 
product space which represents the Hilbert space 

of S combined with an n-level system. We can define a 


map $ 0 /„ operating on the combined system by linear 
extension of the relation ($ 0 dn)(A 0 B) = ($A) 0 B. 
This map <i> 0 /„ thus describes a quantum operation of 
the composite system which acts nontrivially only on the 
first factor representing subsystem S. One then defines 
the map 4> to be n-positive if $ 0 /„ is a positive map, 
and completely positive if $ 0 J„ is a positive map for all 
n. We note that the existence of maps which are positive 
but not completely positive is closely connected to the 
existence of entangled states. We note further that posi¬ 
tivity is equivalent to 1-positivity, and that for a Hilbert 
space with finite dimension Ns = dim'Hs complete pos¬ 
itivity is equivalent to A^ 5 -positivity. 

If the time parameter t now varies over some time in¬ 
terval from 0 to T, where T may be finite or infinite, we 
obtain a one-parameter family of dynamical maps. 


$ = {$* |0<ts^r,$o = /}, (10) 

where I denotes the unit map, and the initial environ¬ 
mental state /9 e( 0) used to construct <I)t is still kept fixed. 
This family contains the complete information on the dy¬ 
namics of all initial open system states over the time in¬ 
terval [0,T] we are interested in. 


2. Divisibility and time-local master equations 

Let us suppose that the inverse of exists for all times 
t ^ 0. We can then define a two-parameter family of 
maps by means of 

^t,s = ^tK\ (11) 

such that we have o = ‘I’t and 

^t.o = $t,s$s.o- (12) 

The existence of the inverse for all positive times thus al¬ 
lows us to introduce the very notion of divisibility. While 
$4 0 and ^>5 0 are completely positive by construction, the 
map need not be completely positive and not even 
positive since the inverse of a completely positive 
map $s need not be positive. The family of dynamical 
maps is said to be P-divisible if is positive, and CP- 
divisible if is completely positive for all t ^ s > 0. 
A simple example for a CP-divisible process is provided 
by a semigroup $t = with a generator C in Lindblad 
form. In this case we obviously have $t_s = which 

is trivially completely positive. As we will discuss later 
on there are many physically relevant models which give 
rise to dynamical maps which are neither CP-divisible 
nor P-divisible. 

An interesting property of the class of processes for 
which exists is given by the fact that they always 
lead to a time-local quantum master equation for the 
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open system states with the general structure 


J^PS=^tPS 


(13) 


= -* [Hsit),ps] 


A,{t)psAl{t) - -{Al{t)A,(t),ps} 


The form of this master equation is very similar to that 
of a Lindblad master equation, where however the Hamil¬ 
tonian contribution Hs{t), the Lindblad operators Ai{t) 
(which can be supposed to be linearly independent), and 
the decay rates 7 i(t) may depend on time since the pro¬ 
cess does not represent a semigroup, in general. Note 


that Eq. (13) involves a time-dependent generator /Ct, 


but no convolution of the open system states with a mem¬ 


ory kernel as in the Nakajima-Zwanzig equation (Naka- 


jima 1958 Zwanzig 1960). Master equations of the form 
(13) can be derived employing the time-convolutionless 


projection operator technique (Chaturvedi and Shibata 


1979 Shibata et al. 1977). 


It is an important open problem to formulate general 
necessary and sufficient conditions under which the mas¬ 


ter equation (13) leads to a completely positive dynamics. 


If the process represents a semigroup, the Hamiltonian 
Hg, the operators Ai and the rates must be time- 
independent and a necessary and sufficient condition for 
complete positivity of the dynamics is simply that all de¬ 
cay rates are positive, ^ 0. This is the famous Gorini- 


Kossakowski-Sudarshan-Lindblad theorem (Gorini et al. 
1976[ [Lindblad 1976). However, we will see later by sev¬ 


eral examples that in the time-dependent case the rates 
7 i(t) can indeed become temporarily negative without 
violating the complete positivity of the dynamics. 

On the other hand, for divisible quantum processes it 
is indeed possible to formulate necessary and sufficient 


conditions. In fact, the master equation (13) leads to a 


CP-divisible dynamics if and only if all rates are positive 
for all times, 7 i(t) ^ 0, which follows from a straight¬ 
forward extension of the Gorini-Kossakowski-Sudarshan- 


Lindblad theorem. Moreover, the master equation (13) 


leads to a P-divisible dynamics if and only if the weaker 
conditions 




(14) 


ing values in a discrete set {cCijigR is characterized 
by a hierarchy of joint probability distributions Pn = 
Pn{xn, tn'j Xn-iAn-i', ■ ■ - IXiAi) for aU n G N and times 
tn ^ tn-i ^ ^ ^ 0. The distribution P„ yields the 

probability that the process takes on the value xi at time 
ti, the value X 2 at time ^ 2 , ■ • ■ j and the value at time 
tn. In order for such a hierarchy to represent a stochastic 
process the Kolmogorov consistency conditions must be 
satisfied which, apart from conditions of positivity and 
normalization, require in particular the relation 

^ ^ (.Xn 1 f n 5 ■ • ■ 5 7 l"m ; • ■ • ; 5 ) 

Xm 

— Pn—l {.Xn (^^) 

connecting the n-point probability distribution Pn to the 
(n — l)-point probability distribution Pn-i- 
A stochastic process X{t) is said to be Markovian if 
the conditional probabilities defined by 

; ^n+1 \ ^n ; ; • ■ • ; ) 

_ -fn+1 (^n+15 ^n+15 ■ • ■ 5 ) 

-fn (^n 1 i ■ • ■ 7^17^1) 

satisfy the relation 

(^n+l itn+l \ Xn^ tn: ■ • ■ : XiAi'j 
— -f^l| 1 (^n-t-15 fn-t-1 fn)- 

This is the classical Markov condition which means that 
the probability for the stochastic process to take the value 
Xn+i at time tn+i, under the condition that it assumed 
values Xi at previous times ti, only depends on the last 
previous value Xn at time In this sense the pro¬ 
cess is said to have no memory, since the past history 
prior to is irrelevant to determine the future once we 
know the value Xn the process assumed at time t„. Note 
that Eq. 0 imposes an infinite number of conditions 
for all n-point probability distributions which cannot be 
checked if only a few low-order distributions are known. 

The Markov condition 0 substantially simplifies the 
mathematical description of stochastic processes. In fact, 
one can show that under this condition the whole hi¬ 
erarchy of joint probability distributions can be recon¬ 
structed from the initial 1-point distribution Pi(a:o,0) 
and the conditional transition probability 


(16) 

(17) 


hold for all orthonormal bases {|n)} of the open system 
and all n ^ m. This statement can be obtained by apply¬ 
ing a characterization of the generators of positive semi¬ 
groups due to Kossakowski (1972a|b). 


B. Classical versus quantum non-Markovianity 

1. Classical stochastic processes and the Markov condition 


T{x,t\y,s) = Pi|i(a;,t|?/,s) 
by means of the relations 

7 tn, Xn—1 7 tn— 1 ^ • ■ • ^ 

n — 1 

= Y]_T{Xi+l,ti+l\Xi:U)Pl{xi,tl) 

and 


(18) 


(19) 


In classical probability theory (Gardiner, 1985 


van 


Kampen 1992) a stochastic process X{t), t ^ 0, tak- 


Pi{xi,ti) = '^T{xi,ti\xo,0)Pi{xo,0). (20) 
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For a Markov process the transition probability has to 
obey the Chapman-Kolmogorov equation 

T{x,t\y,s) = ^T{x,t\z,T)T{z,T\y,s) (21) 


for t ^ T ^ s. Thus, a classical Markov process is 
uniquely characterized by a probability distribution for 
the initial states of the process and a conditional tran¬ 
sition probability satisfying the Chapman-Kolmogorov 
equation (211. Indeed, the latter provides the necessary 
condition in order to ensure that the joint probabilities 
defined by (19) satisfy the condition (15), so that they 
actually define a classical Markov process. Provided the 
conditional transition probability is differentiable with 
respect to time (which will always be assumed here) one 
obtains an equivalent differential equation, namely the 
Chapman-Kolmogorov equation: 


jT{x,t\y,s) (22) 

= X! \Wxzit)T{z, t\y, s) - W:,x{t)T{x, t\y, s) , 

Z 


where Wzx{t) ^ 0, representing the rate (probability per 
unit of time) for a transition to the state z given that 
the state is x at time t. An equation of the same struc¬ 
ture holds for the 1-point probability distribution of the 
process: 



{x,t) = Y, [Wxzit)Pi{z,t)-Wzxit)Piix,t)\, (23) 


which is known as Pauli master equation for a classical 
Markov process. The conditional transition probability 
of the process can be obtained solving Eq. (23) with the 
initial condition Pi{x, s) = S^y 


2. Non-Markovianity in the quantum regime 


The above definition of a classical Markov process can¬ 
not be transferred immediately to the quantum regime 


(Vacchini et al. 2011). In order to illustrate the arising 


difficulties let us consider an open quantum system as 
described in Sec. |II.A[ Suppose we carry out projective 
measurements of an open system observable X at times 
tn ^ tn -1 ^ ^ ^ 0. For simplicity we assume that 

this observable is non-degenerate with ^ectral decom¬ 
position X = write the 

unitary evolution superoperator as UtpsE = UtPssUt 
and the quantum operation corresponding to the mea¬ 
surement outcome x as MxPse = \p>x){p>x\psE\yyx){^x\- 
Applying the Born rule and the projection postulate one 
can then write a joint probability distribution 

Pn {Xn , i-n : Xji—i , t^i— l,...,:ri,ti) 

= tr {Mx„Ut^-t^_, ... Mx,Ut,psE{0)} (24) 


which yields the probability of observing a certain se¬ 
quence Xn, Xn- 1 ,. ■ ■ ,xi of measurement outcomes at the 
respective times tn, tn-i ,..., ti. Thus, it is of course pos¬ 
sible in quantum mechanics to define a joint probabil¬ 
ity distribution for any sequence of measurement results. 
However, as is well known these distribution do in gen¬ 


eral not satisfy condition (15) since measurements carried 


out at intermediate times in general destroy quantum in¬ 
terferences. The fact that the joint probability distribu¬ 


tions (24) violate in general the Kolmogorov condition 


(15) is even true for closed quantum system. Thus, the 


quantum joint probability distributions given by (24) do 


not represent a classical hierarchy of joint probabilities 
satisfying the Kolmogorov consistency conditions. More 
generally, one can consider other joint probability distri¬ 
butions corresponding to different quantum operations 
describing non-projective, generalized measurements. 

For an open system coupled to some environment mea¬ 
surements performed on the open system not only influ¬ 
ence quantum interferences but also system-environment 
correlations. For example, if the system-environment 
state prior to a projective measurement at time ti is given 
by psE{ti), the state after the measurement conditioned 
on the outcome x is given by 


p'sEi^i) — 


MxPSEjtj) 

trMxPSEiU) 


= \Px){p>x\‘^ P%{ti), (25) 


where p% is an environmental state which may depend 
on the measurement result x. Hence, projective mea¬ 
surements completely destroy system-environment corre¬ 
lations, leading to an uncorrelated tensor product state 
of the total system, and, therefore, strongly influence the 
subsequent dynamics. 

We conclude that an intrinsic characterization and 
quantification of memory effects in the dynamics of open 
quantum systems, which is independent of any prescribed 
measurement scheme influencing the time evolution, has 
to be based solely on the properties of the dynamics of 
the open system’s density matrix ps{t). 


C. Quantum non-Markovianity and information flow 


The first approach to quantum non-Markovianity to 
be discussed here is based on the idea that memory ef¬ 
fects in the dynamics of open systems are linked to the 
exchange of information between the open system and 
its environment: While in a Markovian process the open 
system continuously looses information to the environ¬ 
ment, a non-Markovian process is characterized by a flow 
of information from the environment back into the open 


system (Breuer et aZ. 2009 Laine et aZ., 2010a). In such a 


way quantum non-Markovianity is associated to a notion 
of quantum memory, namely information which has been 
transferred to the environment, in the form of system- 
environment correlations or changes in the environmen- 
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tal states, and is later retrieved by the system. To make 
this idea more precise we employ an appropriate distance 
measure for quantum states. 


1. Trace distance and distinguishability of quantum states 


The trace norm of a trace class operator A is defined by 
||A|| = tr|yl|, where the modulus of the operator is given 
by \A\ = V A^ A. If A is selfadjoint the trace norm can be 
expressed as the sum of the moduli of the eigenvalues 
of A counting multiplicities, ||A|| = |ai|. This norm 

leads to a natural measure for the distance between two 
quantum states and p^ known as trace distance: 

D{p\p^) = \\\p^-p% (26) 


The trace distance represents a metric on the state space 
S{'H) of the underlying Hilbert space T-L. We have the 
bounds 0 ^ D{p^,p^) < 1, where D{p^,p^) = 0 if and 
only if and D{p^,p'^) = 1 if and only if p^ and 

p^ are orthogonal. Recall that two density matrices are 
said to be orthogonal if their supports, i.e. the subspaces 
spanned by their eigenstates with nonzero eigenvalue, are 
orthogonal. The trace distance has a series of interest¬ 
ing mathematical and physical features which make it a 
useful distance measure in quantum information theory 


(Nielsen and Chuang 20001. There are two properties 


which are the most relevant for our purposes. 

The first property is that the trace distance between 
two quantum states admits a clear physical interpreta¬ 
tion in terms of the distinguishability of these states. To 
explain this interpretation consider two parties, Alice and 
Bob, and suppose that Alice prepares a quantum system 
in one of two states, p^ or p^, with a probability of | 
each, and sends the system to Bob. The task of Bob 
is to find out by means of a single quantum measure¬ 
ment whether the system is in the state p^ or p^. It can 
be shown that the maximal success probability Bob can 
achieve through an optimal strategy is directly connected 
to the trace distance: 


Rmax = ^ [1 -b£>(p\p^)] . 


(27) 


Thus, we see that the trace distance represents the bias 
in favor of a correct state discrimination by Bob. For this 
reason the trace distance D{p^,p^) can be interpreted as 
the distinguishability of the quantum states p^ and p^. 
For example, suppose the states prepared by Alice are 
orthogonal such that we have D{p^,p^) = 1. In this 
case we get Pmax = 1) which is a well known fact since 
orthogonal states can be distinguished with certainty by 
a single measurement, an optimal strategy of Bob being 
to measure the projection onto the support of p^ or p^. 

The second important property of the trace distance is 
given by the fact that any completely positive and trace 



FIG. 2 (Color online) Illustration of the loss of distinguisha¬ 
bility quantified by the trace distance D between density ma¬ 
trices as a consequence of the action of a quantum dynamical 
map <l>t. Alice prepares two distinct states, but a dynamical 
map acts as a noisy channel, thus generally reducing the in¬ 
formation available to Bob in order to distinguish among the 
states by performing measurements on the system only. 


preserving map A is a contraction for the trace distance, 
i.e. we have 


D{Ap\Ap^)^D{p\p^) 


(28) 


for all states p^’^. In view of the interpretation of the 
trace distance we thus conclude that a trace preserving 
quantum operation can never increase the distinguisha¬ 
bility of any two quantum states. We remark that the 
equality sign in (28) holds if A is a unitary transforma¬ 
tion, and that (28) is also valid for trace preserving maps 
which are positive but not completely positive. 


2. Definition and quantification of memory effects 


In the preceding subsection we have interpreted the 
trace distance of two quantum states as the distinguisha¬ 
bility of these states, where it has been assumed that 
the quantum state Bob receives is identical to the state 
prepared by Alice. Suppose now that Alice prepares 
her states P 5 ^( 0 ) as initial states of an open quantum 
system S coupled to some environment E. Bob will 
then receive at time t the system in one of the states 
where denotes the corresponding 
quantum dynamical map. This construction is equiva¬ 
lent to Alice sending her states through a noisy quantum 
channel described by the completely positive and trace 


preserving map . According to (28) the dynamics gen¬ 


erally diminishes the trace distance and, therefore, the 
distinguishability of the states. 


^(Ps(0:Ps(0) ^ ^(Ps(0): 7's'(0))i 


(29) 


such that it will in general be harder for Bob to discrimi¬ 
nate the states prepared by Alice. This fact is illustrated 
schematically in Fig.|^ in which the distinguishability of 
states prepared by Alice decreases due to the action of a 
dynamical map, so that Bob is less able to discriminate 
among the two. Thus, we can interpret any decrease of 
the trace distance D{pg{t), Ps{t)) as a loss of information 


















from the open system into the environment. Conversely, 
if the trace distance D{p\{t), pg{t)) increases, we say that 
information flows back from the environment into the 
open system. 

This interpretation naturally leads to the following def¬ 
inition: A quantum process given by a family of quantum 
dynamical maps is said to be Markovian if the trace 
distance D{pg{t), pg{t)) corresponding to all pair of ini¬ 
tial states Pg(0) and p|(0) decreases monotonically for all 
times t ^ 0. Quantum Markovian behavior thus means 
a continuous loss of information from the open system 
to the environment. Conversely, a quantum process is 
non-Markovian if there is an initial pair of states p;s(0) 
and p|(0) such that the trace distance D{pg{t), p'g{t)) 
is non-monotonic, i.e. starts to increase for some time 
t > 0. If this happens, information flows from the envi¬ 
ronment back to the open system, which clearly expresses 
the presence of memory effects: Information contained 
in the open system is temporarily stored in the environ¬ 
ment and comes back at a later time to influence the 
system. A crucial feature of this definition of quantum 
Markovian process is the fact that non-Markovianity can 
be directly experimentally assessed, provided one is able 
to perform tomographic measurement of different initial 
states at different times during the evolution. No prior 
information on the dynamical map d>( is required, apart 
from its very existence. 


the total system states minus the distinguishability of the 
open system states at time t. In other words, Iext{t) can 
be viewed as the gain in the state discrimination Bob 
could achieve if he were able to carry out measurements 
on the total system instead of measurements on the open 
system only. We can therefore interpret Xint(i) as the 
information inside the open system and lext(i) as the 
amount of information which lies outside the open sys¬ 
tem, i.e., as the information which is not accessible when 
only measurements on the open system can be performed. 
Obviously, we have Xint(i) ^ 0 and lext(i) ^ 0. Since the 
total system dynamics is unitary the distinguishability of 
the total system states is constant in time. Moreover we 
have H(p^£;( 0),/9|£;(0)) = £)(p;^(0),p|(0)) because, by 
assumption, the initial total states are uncorrelated with 
the same reduced environmental state. Hence, we obtain 


Xint(t) -klext(0 =Xint(0) = COnst. 


(32) 


Thus, initially there is no information outside the open 
system, Xext(O) = 0. If Xint(t) decreases Xext(t) must 
increase and vice versa, which clearly expresses the idea 
of the exchange of information between the open system 
and the environment illustrated in Fig. Employing 
the properties of the trace distance one can derive the 
following general inequality ( Laine et al.\ |20I0b ) which 
holds for alH ^ 0: 



FIG. 3 (Color online) Illustration of the information flow 
between an open system and its environment according to 
Eq. (321. Left: The open system looses information to the en¬ 
vironment, corresponding to a decrease of Xint(t) and Marko¬ 
vian dynamics. Right: Non-Markovian dynamics is charac¬ 
terized by a backflow of information from the environment to 
the system and a corresponding increase of Xint(t). 


In order to explain in more detail this interpretation 
in terms of an information flow between system and en¬ 
vironment let us define the quantities 

lintit) = D{pUt),pUt)) (30) 

and 

Xext(i) = D{pg^{t), pg^{t)) — D{pg{t), Ps(t))- (31) 

Here, Xint(t) is the distinguishability of the open system 
states at time t, while Xgxt (t) is the distinguishability of 


Xcxt(t)< ^(PsE(i), Ps(i) ^ PE(t)) (33) 

+-D(4ij(t), Ps(t) O p%(t)) + n(plj(t), pl(t)). 

The right-hand side of this inequality consists of three 
terms: The first two terms provide a measure for the cor¬ 
relations in the total system states pg^(t), given by the 
trace distance between these states and the product of 
their marginals, while the third term is the trace distance 
between the corresponding environmental states. Thus, 
when Tcxt(t) increases over the initial value Xoxt(O) = 0 
system-environment correlations are built up or the envi¬ 
ronmental states become different, implying an increase 
of the distinguishability of the environmental states. This 
clearly demonstrates that the corresponding decrease of 
the distinguishability Xint(t) of the open system states 
always has an impact on degrees of freedom which are 
inaccessible by measurements on the open system. Con¬ 
versely, ifXint(t) starts to increase at some point of time t, 
the corresponding decrease of Xgxt (i) implies that system- 
environment correlations must be present already at time 
t, or that the environmental states are different at this 
point of time. 

On the ground of the above definition for non- 
Markovian dynamics one is naturally led to introduce 
the following measure for the degree of memory effects 

A/’($) = max f dt a{t), (34) 

Ja>0 
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where 




(35) 


denotes the time derivative of the trace distance of the 


evolved pair of states. In Eq. (34) the time integral is 


extended over all intervals in which a{t) > 0, i.e., in 
which the trace distance increases with time, and the 
maximum is taken over all pairs of initial states of 
the open system’s state space S{'Hs)- Thus, A/'(<i)) is a 
positive functional of the family of dynamical maps 
which represents a measure for the maximal total flow 
of information from the environment back to the open 
system. By construction we have = 0 if and only 

if the process is Markovian. 

The maximization over all pairs of quantum states in 


Eq. (34) can be simplified substantially employing sev¬ 


eral important properties of the functional Af {^) and 
the convex structure of the set of quantum states. A 
certain pair of states pg'^ is said to be an optimal state 
pair if the maximum in Eq. (34) is attained for this pair 


of states. Thus, optimal state pairs lead to the maxi¬ 
mal possible backflow of information during their time 
evolution. One can show that optimal states pairs 
lie on the boundary of the state space, and are in par¬ 


ticular always orthogonal (Wifimann et al.. 2012). This 


is a quite natural result in view of the interpretation in 
terms of an information flow since orthogonality implies 
that D{pg, p|) = 1, which shows that optimal state pairs 
have maximal initial distinguishability, corresponding to 
a maximal amount of initial information. An even more 
drastic simplification is obtained by employing the fol¬ 
lowing equivalent representation (Liu et al. 2014): 


A/’(4’) = max / dt aU), 

P^dU{po)Js^O 


where 


a{t) = 


dt 


D{^tP, ^tPo) 


D ip,Po) 


(36) 


(37) 


is the time derivative of the trace distance at time t di¬ 


vided by the initial trace distance. In Eq. (36) po is a 


fixed point of the interior of the state space and dU{po) 
an arbitrary surface in the state space enclosing, but not 
containing po. The maximization is then taken over all 
points of such an enclosing surface. This representation 
often significantly simplifies the analytical, numerical or 
experimental determination of the measure since it only 
involves a maximization over a single input state p. It 
is particularly advantageous if the open system dynam¬ 
ics has an invariant state and if this state is taken to be 
Po, such that only one state of the pair evolves in time. 


Equation (36) may be called local representation for it 


shows that optimal state pairs can be found in any lo¬ 
cal neighborhood of any interior point of the state space. 


Eurthermore, the representation reveals the universality 
of memory effects, namely that for a non-Markovian dy¬ 
namics memory effects can be observed everywhere in 
state space. 


3. Generalizing the trace distance measure 

There is an interesting generalization of the above defi¬ 
nition and quantification of non-Markovian quantum dy¬ 


namics first suggested by Chruscihski et al. (2011). Re¬ 


turning to the interpretation of the trace distance de¬ 


scribed in Sec. II.C.l we may suppose that Alice pre¬ 
pares her quantum system in the states pg or pg with 
corresponding probabilities pi and P 2 , where P 1 +P 2 = 1- 
Thus, we assume that Alice gives certain weights to her 
states which need not be equal. In this case one can again 
derive an expression for the maximal probability for Bob 
to identify the state prepared by Alice: 


-Pmax — X (1 + WPiPs ~ P2Ps\ 


(38) 


As follows from (38) the quantity ||piPg — P 2 Ps\\ gives 
the bias in favor of the correct identification of the state 
prepared by Alice. We note that th e operator A = p ipg — 
P2p‘s known as Helstrom matrix (Helstrom 1967), and 
that Eq. ([38|) reduces to Eq. (p^ in the unbiased case 


Pi = P 2 = Note that the interpretation in terms 
of an information flow between system and environment 
still holds for the biased trace distance, as can be seen 


replacing in (30) and (31) the trace distance by the norm 


of the Helstrom matrix. 

Following this approach we can define a process to be 
Markovian if the function ||$t(pipg — P 2 p|)|| decreases 
monotonically with time for all pi _2 and all pairs of ini¬ 
tial open system states p^^. Consequently, the extended 
measure for non-Markovianity is then given by 


where 


AC)®) = max / dt cr(t). 
Pi ,P% ^ (7>0 


= ^W^tiPlPs - P2 pI)\\- 


(39) 


(40) 


This generalized definition of quantum non-Markovianity 
leads to several important conclusions. First, we note 
that also for the measure ( [3^ optimal state pairs are or¬ 
thogonal, satisfying HpiPg — P 2 Psl| = 1 and, therefore, 
corresponding to maximal initial distinguishability. The 


maximum in (39) can therefore be taken over all Hermi- 


tian (Helstrom) matrices A with unit trace norm. This 


leads to a local representation analogous to Eq. (36): 


A7($) = max [ dt a{t), (41) 

AedUio) 
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where a{t) = ^||$tA||/||A|j, and dU{0) is a closed sur¬ 
face which encloses the point A = 0 in the K-linear vector 
space of Hermitian matrices. 



FIG. 4 (Color online) Schematic picture of the relations be¬ 
tween the concepts of quantum Markovianity and divisibility 
of the family of quantum dynamical maps The blue area 
to the left of the thick black line represents the set of Marko¬ 
vian quantum processes for which the inverse of the dynam¬ 
ical map exists, which is identical to the set of P-divisible 
processes. This set contains as subset the set of CP-divisible 
processes, depicted in green. The blue area to the right repre¬ 
sents the non-Markovian processes which are therefore neither 
CP nor P-divisible. The zig-zag line dividing the two regions 
points to the fact that neither of them is convex. The red 
areas mark the Markovian or non-Markovian processes for 
which the inverse of "Ft does not exist. 


Another important co nclusion is obtained if we as¬ 
sume, as was done in Sec. II.A.2 that exists. It can 
be shown that under this condition the quantum process 
is Markovian if and only if is P-divisible, which follows 
immediately from a theorem by Kossakowski| ( |1972a|b I. 
The relationships between Markovianity and divisibility 
of the dynamical map are illustrated in Fig.|^ where also 
the situation in which the inverse of as a linear map 
does not exist for all times has been considered. The 
generalized characterization of non-Markovianity is sen¬ 
sitive to features of the open system dynamics neglected 
by the trace distance measure (Liu et al. 2013b). It is 
worth stressing that in keeping with the spirit of the ap¬ 
proach, this generalized measure of non-Markovianity, in 
addition to its clear connection with the mathematical 
property of P-divisibility of the dynamical map, can still 
be directly evaluated by means of experiments. More 
specifically, both the quantification of non-Markovianity 
according to the measure (34), as well as its generaliza¬ 
tion (39) can be obtained from the very same data. 


4. Connection between quantum and classical non-Markovianity 


The definition of quantum non-Markovianity of 


Sec. II.C.3 allows to establish an immediate general 
connection to the classical definition. Suppose as in 
Sec. |II.A.2| that the inverse of the dynamical maps 
exists such that we have the time-local quantum mas¬ 


ter equation (13), in which for the sake of simplicity we 


consider only the rates to be time dependent. Suppose 
further that maps operators which are diagonal in a 
fixed basis {|n)} to operators diagonal in the same ba¬ 
sis. If ps(0) is an initial state diagonal in this basis, then 
psit) can be expressed at any time in the form: 


Psit) = ^P„(t)|n)(n|, 


(42) 


where P„(t) denote the time-dependent eigenvalues. The 
quantum master equation then leads to the following 
equation of motion for the probabilities Pn{t): 

[Wnmit)Pmit) - (t)P„ (t)] , (43) 


where 


Wnm{t) = y^7i(t)|(n|Ai|m) 


(44) 


We observe that Eq. (43) has exactly the structure of a 


classical Pauli master equation (23). Therefore, Eq. (43) 


can be interpreted as differential Chapman-Kolmogorov 


equation (22) for the conditional transition probability of 


a classical Markov process if and only if the rates (44) are 


positive. As we have seen a quantum process is Marko¬ 
vian if and only if the dynamics is P-divisible. Accord¬ 


ing to condition (14) P-divisibility implies that the rates 


(44) are indeed positive. Thus, we conclude that any 


quantum Markovian process which preserves the diago¬ 
nal structure of quantum states in a fixed basis leads to a 
classical Markovian jump process describing transitions 
between the eigenstates of the density matrix. 


D. Alternative approaches 

In Sec. |II.C| we have considered an approach to quan¬ 
tum non-Markovianity based on the study of the dynam¬ 
ical behavior of the distinguishability of states prepared 
by two distinct parties, Alice and Bob. This viewpoint 
directly addresses the issue of the experimental detection 
of quantum non Markovian processes, which can be ob¬ 
tained by suitable quantum tomographic measurements. 
It further highlights a strong connection between quan¬ 
tum non-Markovianity and quantum information theory, 
or more specifically quantum information processing. 

Alternative approaches to the study of quantum non- 
Markovian dynamics of open systems have also been re¬ 
cently proposed. We shall try to briefly discuss those 
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which are more closely related to the previous study, 
without any claim to properly represent the vast liter¬ 
ature on the subject. For a review on recent results 


on quantum non-Markovianity see Rivas et al. (20141. 


The different approaches can to some extent be grouped 
according to whether they propose different divisibility 


properties of the quantum dynamical map (Chruscihski 


and Maniscalcol 2014| Hall et al.[ 20141|Hou et al. 

2011 

Rivas et al. 2010 Wolf et al. 2008), different quantifiers 

of the distinguishability between states (Chruscihski and 

Kossakowskil 2012| Dajka et al.. 2011 |Lu et o/.| [20101 

Vasile et al. 

2011a Wifimann et aZ. 2013) or the st 

udy of 
nation 
otonic 

other quant 
concepts or 
and an osci 

ities, be they related to quantum inforn 
not, which might exhibit both a mon 

lating behavior in time (Bylicka et al. 

2014 

Fanchini et al.\ 2014| [Haseli et a/.[[2014| Lu et al. 

2010 

Luo et al. 

m2\. 


Pl{t) and Pi{t). Therefore these two equivalent prop¬ 
erties capture the feature of a classical Markovian pro¬ 
cess at the level of the one-point probability distribution, 
which is all we are looking for in order to obtain an in¬ 
trinsic definition of quantum non-Markovianity. 

Relying on the fact that in the theory of open quantum 
systems complete positivity suggests itself as a natural 
counterpart of positivity, as an alternative approach to 
quantum non-Markovianity it has been proposed to iden¬ 
tify quantum Markovian dynamics with those dynamics 
which are actually CP-divisible. This criterion was first 


suggested by Rivas et al. (2010). In order to check CP- 


1. Quantum non-Markovianity and CP-divisibiiity 


Let us first come back to the notion of divisible maps 


divisibility of the time evolution one actually needs to 
know the mathematical expression of the transition maps 
or of the infinitesimal generator K,f Indeed given the 
maps one can determine their complete positivity 
studying the associated Choi matrices, while as discussed 
in Sec. EH CP-divisibility corresponds to positivity at 
all times of the rates appearing in (13). In order to 
quantify non-Markovianity as described by this criterion, 
different related measures have been introduced, relying 
on an estimate of the violation of positivity of the Choi 


introduced in Sec. |II.A.2[ which we express in the form matrix 1 

Rivas et al. 

2010 

or on various quantifiers of 

the negativity of the rates 

Hall et al. 2014 

). In particu- 


- 41^ 7-4^7- 3, t^T ^ s ^ 0, 


(45) lar, Rivas et al. (2010) suggest to consider the following 


where the maps have been defined according to (11). 


measure of non-Markovianity based on CP-divisibility 


In the classical case considering the matrices whose el¬ 
ements are given by the conditional transition prob- 


.A'rhp(4>) = / dtg{t) 


abilities (18) according to {h.t^s)xy = T{x,t\y,s) the 
Chapman-Kolmogorov equation takes the form 


with 




t> T ^ Q, 


(46) 


g{t) = lim 
€—>■0 + 


®J) (IT)(4/1)11-1 


(48) 


(49) 


where each Aj ^ is a stochastic matrix, which obviously 


bears a strict relationship with (45), once one consid¬ 


ers that positive maps are natural quantum counter¬ 
part of classical stochastic matrices. Given the fact that 
in the classical case the Markov condition 0 entails 
the Chapman-Kolmogorov equation, and building on the 
conclusion drawn at the end of Sec. |II.B.2[ i.e. that 
an intrinsic characterization of memory effects has to be 
based on the system’s density operator only, whose dy¬ 
namics is fully described in terms of the time evolution 
maps, one is led to associate quantum Markovianity with 
P-divisibility. This observation reinforces the result ob¬ 
tained in Sec. |II.C.3| relying on the study of the informa¬ 
tion flow as quantified by the generalized trace distance 
measure. Actually it can be shown that the P-divisibility 


where |T) is a maximally entangled state between open 
system and ancilla. 

A further feature which can be considered in this con¬ 
text, and which has been proposed in [Lorenzo et al. 
(2013) as another signature of non-Markovianity is the 


dynamical behavior of the volume of the accessible states, 
according to their parametrization in the Bloch repre¬ 
sentation, identifying non-Markovianity with the growth 
of this volume, which also allows for a geometrical de¬ 
scription. This study can be performed e.g. upon knowl¬ 
edge of the rates appearing in the time-local master equa¬ 
tion, and appear to be a strictly weaker requirement for 


Markovianity with respect to P-divisibility (Chruscihski 


and Wudarski 2013). 


condition (461 is equivalent to the monotonicity property 


of the function 


K{Pl{t),Pi{t)) = \piPl{x,t) - p2Pi{x,t)\ (47) 


which provides a generalization of the Kolmogorov dis¬ 
tance between two classical probability distributions 


To explain the relation between P-divisibility and CP- 
divisibility suppose that Alice prepares states or p| on 
an extended space 'Hs®'Ps and sends the states through 
the channel Tt = TtC/ to Bob. The dynamical map T* is 
then Markovian, in the sense of P-divisibility, if and only 
if is CP-divisible. Note however that this construction 
requires the use of a physically different system, described 
by the tensor product space Hs ® P-S, which carries the 
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information transferred from Alice to Bob. More specifi¬ 
cally, the maximization over the initial states must allow 
for entangled states which clearly explains the differ¬ 
ence between the definitions of non-Markovianity for 
and $t. This implies in particular that an experimental 
assessment of non-Markovianity according to the crite¬ 
rion of CP-divisibility would in general imply the capabil¬ 
ity to prepare entangled states and perform tomographic 
measurements on the extended space. Alternatively one 
should be able to perform process tomography instead of 
state tomography, with an analogous scaling with respect 
to the dimensionality of the system. 


2. Monotonicity of correlations and entropic quantities 


Another line of thought about non-Markovianity of a 
quantum dynamics is based on the study of the behavior 
in time of quantifiers of correlations between the open 
system of interest and an ancilla system. In this respect 


both entanglement (Rivas et al. 2014) and quantum mu¬ 
tual information have been considered (Luo et al. 20121. 


In this setting one considers an ancilla system and stud¬ 
ies the behavior in time of the entanglement or the mu¬ 
tual information of an initially correlated joint system- 
ancilla state. Given the fact that both quantities are 
non-increasing under the local action of a completely pos¬ 
itive trace preserving transformation, CP-divisibility of 
the time evolution would lead to a monotonic decrease 
of these quantities. Their non monotonic behavior can 
therefore be taken as a signature or a definition of quan¬ 


tum non-Markovianity. Actually, Rivas et al. (2014) have 


considered CP-divisibility as the distinguishing feature 
of non-Markovianity, so that a revival of entanglement 
in the course of time is just interpreted as a witness of 


non-Markovianity. In (Luo et al.. 20121 instead failure 


of monotonicity in the loss of quantum mutual informa¬ 
tion is proposed as a new definition of non-Markovianity. 
In both cases one considers as initial state a maximally 
entangled state between system and ancilla. The quan¬ 
tification of the effect is obtained by summing up the in¬ 
creases over time of the considered quantity. Still other 
approaches have connected non-Markovianity with non 
monotonicity of the quantum Fisher information, consid¬ 
ered as a quantifier of the information flowing between 
system and environment. In this case however the infor¬ 
mation flow is not directly traced back to a distance on 


the space of states (Lu et al.. 2010). 


In the trace distance approach one studies the behav¬ 
ior of the statistical distinguishability between states, 
and relying on the introduction of the Helstrom matrix 


(Chruscihski et al. 20111 this notion can be directly re¬ 


lated to a divisibility property of the quantum dynamical 
map. The generalization of the trace distance approach 
also allows for a stronger connection to a quantum infor¬ 
mation theory viewpoint. Indeed one can consider the 


quantum dynamical maps as a collection of time depen¬ 
dent channels. In this respect the natural question is 
how well Bob can recover information on the state pre¬ 
pared by Alice performing a measurement after a given 
time t. A proposal related to this viewpoint has been 
considered in (Bylicka et al. 2014), where two types 
of capacity of a quantum channel, namely the so-called 
entanglement-assisted classical capacity and the quan¬ 
tum capacity, have been considered. While the quantum 
data processing inequality (Nielsen and Chuang 2000) 
warrants monotonicity of these quantities provided the 
time evolution is CP-divisible, lack of this monotonicity 
has been taken as definition of non-Markovianity of the 
dynamics. The connection between non-Markovianity 
and the quantum data processing inequality has also been 
recently studied by Buscemi and Datta (2014). 


III. MODELS AND APPLICATIONS OF 
NON-MARKOVIANITY 

In this section we first address the description of simple 
prototypical systems, for which one can exactly describe 
the non-Markovian features of the dynamics according 
to the concepts and measures introduced in Sec. [ll] In 
particular we will stress how a non-Markovian dynamics 
leads to a recovery of quantum features according to the 
notion of information backflow. In the second part we 
present applications of the introduced non-Markovianity 
measures to a number of more involved models illus¬ 
trating how memory effects allow to analyze their dy¬ 
namics and can be related to characteristic properties of 
complex environments, thus suggesting to use the non- 
Markovianity of the open system dynamics as a quantum 
probe for features of the environment. 


A. Prototypical model systems 


We consider first a two-level system undergoing pure 
decoherence. This situation can be described by a time- 


convolutionless master equation of the form (13) with 


a single decay channel and allows to demonstrate in a 
simple manner how the properties of the environment 
and the system-environment coupling strength influence 
the transition from a Markovian to a non-Markovian dy¬ 
namics for the open quantum system. A similar anal¬ 
ysis is then performed for a dynamics which also takes 
into account dissipative effects, further pointing to the 
phenomenon of the failure of divisibility of the quantum 
dynamical map. We further discuss a dynamics driven 
by different decoherence channels, which allows to bet¬ 
ter discriminate between different approaches to the de¬ 
scription of non-Markovian behavior. Finally, we present 
some results on the non-Markovian dynamics of the spin- 
boson model. 
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1. Pure decoherence model 


Let us start considering a microscopic model for a pure 
decoherence dynamics which is amenable to an exact so¬ 
lution and describes a single qubit interacting with a 
bosonic reservoir. The total Hamiltonian can be writ¬ 
ten as in Eq. (§ with the system Hamiltonian Hs and 
the environmental Hamiltonian He given by 

= ^^00-2, He = '^u}kalak, (50) 

k 

where ujq is the energy difference between ground state 
|0) and excited state |1) of the system, tr^ denotes a Pauli 
matrix, Ofc and al are annihilation and creation opera¬ 
tors for the bosonic reservoir mode labelled by k with 
frequency ujk, obeying the canonical commutation rela¬ 
tions [ak,a\,] = 5kk'■ The interaction term is taken to 
be 

Hi^^a^ {guak + 5 ^ 4 ) ( 51 ) 

k 

with coupling constants gk ■ Considering a factorized ini¬ 
tial condition with the environment in a thermal state at 
inverse temperature 13 the model can be exactly solved 
leading to a quantum dynamical map which leaves 
the populations invariant and modifies the off-diagonal 
matrix elements according to 


Pii(l) — Pii(O)) Pio(l) — G(t)pio(0), (52) 


where Pij{t) = {i\ps{t)\j) denote the elements of the in¬ 
teraction picture density matrix psit)- The function G{t) 
is often called decoherence function and in the present 
case it is real since we are working in the interaction pic¬ 
ture. It can be expressed in the form 


G{t) = exp 



du! J{ui) coth 



1 —cos {(jjt) 


, (53) 


where we have introduced the spectral density J{lo) = 
J2k — Wfc) which keeps track of the features of 

the environment relevant for the reduced system descrip¬ 
tion, containing informations both on the environmental 
density of the modes and on how strongly the system 
couples to each mode. It can be shown that the inter¬ 
action picture operator ps{t) obeys a master equation of 


the form (131 with a single Lindblad operator. 


^Ps(i) = lit) [a^ps{t)a:, - psit )], (54) 


where the time dependent decay rate 'y{t) is related to 
the decoherence function G(t) by 




Note that for a generic microscopic or phenomenological 
decoherence model the decoherence function is in general 
a complex quantity even in the interaction picture, and 
in this case the modulus of the function should be consid¬ 


ered in formula (551. The map <I>t is completely positive 


since G(t) ^ 1 which is equivalent to the positivity of the 
time integral of the decay rate r{t) = dt' In or¬ 
der to discuss the non-Markovianity of the obtained time 
evolution as discussed in Sec. III. Cl and Sec. III. D. II we 
consider the divisibility property of the time evolution. 


The master equation (54) has a single channel with decay 


rate 7 (t) and thanks to the strict positivity of the deco¬ 
herence function always exists. However, according 
to the criterion given by Eq. (14), P-divisibility fails at 


times when the decoherence rate becomes negative, which 
is exactly when CP-divisibility is lost. This is generally 
true when one has a master equation in the time-local 


form (13) with a single Lindblad operator, since in this 


case the condition for CP-divisibility, given by the re¬ 
quirement of positivity of the decay rate 7 (t), coincides 


with the condition (14) for P-divisibility. 


The trace distance between time evolved states corre¬ 
sponding to the initial conditions is given according 
to (1^ by 


D <^^pI) = v'a2+G2(t)|6|2, 


(56) 


where a = p\i—pii and b = Pio~Pio denote the difference 
of the populations and of the coherences of the initial 
states respectively. Its time derivative corresponding to 


the quantity (35) in terms of which the non-Markovianity 
measure can be constructed then reads 


rit) = 


G(t)|6|^ 


Va2 + G2(t)|6|2dt 


Git). 


(57) 


The non-monotonicity in the behavior of the trace dis¬ 
tance which is used as a criterion for non-Markovianity 
in Sec. |H.C.2| is therefore determined by the non mono¬ 
tonic behavior of the decoherence function. In order to 


evaluate the measure Eq. (34) one has to maximize over 


all pairs of initial states. As shown by [Wifimann et al. 
(2012) and discussed in Sec. H.C.2 the latter can be re¬ 


stricted to orthogonal state pairs. For the case at hand it 


follows from Eq. (57) that the maximum is obtained for 


antipodal points on the equator of the Bloch sphere, so 
that a = 0 and |6| = 1. The non-Markovianity measure 
(34) is then given by 


A7(4>) = ^[G(4)-G(f4 


(58) 


where t^ and t{. denote initial and final time point of 
the A:-th interval in which G(t) increases. Since the in¬ 
crease of G(t) coincides with the negativity of the de¬ 
coherence rate according to Eq. ( |5^ , the growth of the 
trace distance is in this case equivalent to breaking both 
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P-divisibility and CP-divisibility. Moreover, the gener¬ 
alized trace distance measure considered in Sec. rrmai 


leads to the same expression (58) for non-Markovianity. 


The measure for non-Markovianity of Eq. (48 1 based 


on CP-divisibility takes in this case the simple expression 


A/RHp(‘k) = / dt 7 (t), 


(59) 


'7<0 


which is the integral of the decoherence rate in the time 
intervals in which it becomes negative. An example of 
an experimental realization of a pure decoherence model 
is considered in Sec. |IV.CT| where the decoherence func¬ 
tion is determined by the interaction between the polar¬ 
ization and the frequency degrees of freedom of a photon. 
In this case the decoherence function is, in general, com¬ 
plex valued, so that the previous formulae hold with G{t) 
replaced by 10(^)1. The geometric measure discussed in 
Sec. |II.DTH which connects non-Markovian behavior with 
the growth in time of the volume of accessible states, 
also leads to the same signature for non-Markovianity: 
The volume can again be identified with the decoherence 
function G{t) and therefore as soon as G{t) increases the 
volume grows. 

As we have discussed, in a simple decoherence model 
non-Markovianity can be identified with the revival of 
coherences of the open system, corresponding to a back- 
flow of information from the environment to the system. 


In this situation the time-local master equation (54) de¬ 


scribing the dynamics exhibits just a single channel, so 
that all discussed criteria of non-Markovianity do coin¬ 
cide. 


2. Two-level system in a dissipative environment 


where the decoherence function is a complex function de¬ 
termined by the spectral density J{uj) of the model. In 
particular, denoting by f(t — ti) the two-point correlation 
function of the environment corresponding to the Fourier 
transform of the spectral density, the function G{t) is de¬ 
termined as the solution of the integral equation 


— G{t) — —J dti f{t — ti)G{ti) 


(62) 


with initial condition C?(0) = 1 . Also in this case it is 
possible to write down the exact master equation obeyed 
by the density matrix ps{t) of the system in the form of 
Eq. (131, which in the interaction picture reads (Breuer 


et al. 19991 


=- ^S{t) [a^,ps] 


(63) 


■lit) 


a-psit)a+ - - {a+a-,psit)} 


where the time dependent Lamb shift is given by S{t) = 
—2Q{G{t)/G{t)), and the decay rate can be written as 
7 (t) = —25i(G'(t)/G'(<)), or equivalently 


lit) = - 


\G{t)\ dt 


|G(t)|. 


(64) 


The analysis of non-Markovianity in the present model 
closely follows the discussion in m .A.l| due to the cru¬ 
cial fact that the master equation (63) still has a single 


Lindblad operator. The expression of the trace distance 
is determined by the decoherence function and, with the 
same notation as in (15^, reads 


D <^tpl) = |G(t)|x/|G(t)|2a2 + |5|2, (65) 

so that one has the time derivative 


We now consider an example of dissipative dynamics 
in which the system-environment interaction influences 
both coherences and populations of a two-state system. 
The environment is still taken to be a bosonic bath, so 
that the free contributions to the Hamiltonian are again 


given by Eq. (50), but one considers an interaction term 


in rotating wave approximation given by 


Hi = '^\ gkcr+ak 




(60) 


where cr_ = | 0 )( 1 | and (T+ = | 1 )( 0 | are the lowering and 
raising operators of the system. Such an interaction de¬ 
scribes, e.g., a two-level atom in a lossy cavity. For a 
factorized initial state with the environment in the vac¬ 
uum state, this model is again exactly solvable, thanks 
to the conservation of the number of excitations. The 
quantum dynamical map transforms populations and 
coherences according to 

Piiit) = |G(t)ppii(0), Pioit) = G(t)pio(0), (61) 


cr(t) = 


2|G(t)pa^ -I- | 6 p d 
V^|G(t)|2a2 + |^^ 


( 66 ) 


Complete positivity of the map is ensured by positivity of 
the integrated decay rate T{t), while the trace distance 
shows a non monotonic behavior when the modulus of 
the decoherence function grows with time. 

A typical expression for the spectral density is given 
by a Lorentzian 


T(w) — 7oA^/27r [(wq + A — w)^ -|- A^] (67) 


of width A and centered at a frequency detuned from the 
atomic frequency ujq by an amount A, while the rate 70 
quantifies the strength of the system-environment cou¬ 
pling. Let us first discuss the case when the qubit is in 
resonance with the central frequency of the spectral den¬ 
sity, so that one has a vanishing detuning, A = 0. The 
decoherence function then takes the form 


G{t) = 


cosh 



-I- — smh 
d 



( 68 ) 
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where d = y/X^ — 270 A. For weak couplings, correspond¬ 
ing to 7 o < A/2, the decoherence function G{t) is a real, 


expression for the decoherence function. For this case 


monotonically decreasing function so that (661 is always 


negative and the dynamics is Markovian. In this case, as 
discussed in Sec. |II.A.2[ the map is invertible and since 
the time-convolutionless form of the master equation has 
just a single channel, it is at the same time P-divisible and 
CP-divisible. Note that in the limit 70 ^ A/2 the rate 
7 (t) becomes time-independent, 7 = 70 , such that the 


master equation (63) is of Lindblad form and leads to a 


Markovian semigroup. However, for larger values of the 
coupling between system and environment, namely for 
7 o > A/2, the function |G(t)| displays an oscillatory be¬ 
havior leading to a non-monotonic time evolution of the 
trace distance and, hence, to non-Markovian dynamics. 
Thus, the trace distance measure of Eq. ( [M| is zero for 
7 o < A/2 and starts at the threshold 70 = A/2 to increase 
continuously to positive values (Laine et al ., 2010a |. It is 


the derivative of the trace distance (661 takes the simple 
form a{t) = — 7 (t) exp[—r(t)] with T{t) = /g 7 (t') dt', 
which puts into evidence the direct connection between 
the direction of information flow and the sign of the decay 
rate. 


3. Single qubit with multiple decoherence channels 


In the preceding examples we have shown that for a 


master equation in the time-convolutionless form (13) 
with a single Lindblad operator A{t), all the different 
considered criteria for non-Markovianity coincide. It is 
therefore of interest to discuss an example in which one 
has multiple decoherence channels with generally differ¬ 
ent rates. As we shall shortly see, such model allows 
to discriminate among the different definitions of non- 
Markovian dynamics. To this end, we take a phenomeno- 


breaks down. We also remark that for this model opti¬ 
mal state pairs still correspond to antipodal points of the 


actly coincides with the point where the perturbation ex- 

logical approach and following Chruscihski and Manis- 

pansion of the time-convolutionless master equation (63) 

calco| 

2014D; 

Chruscihski and Wudarski 

(2013 

1 ; and 

Vac- 


chini ( 2012 ) we consider for a two-level system the master 


equation 


equator of the Bloch sphere (Xu et al. 20101, and that 


the criterion for non-Markovianity based on the Helstrom 
matrix considered in Sec. III. C. P leads to the same results. 

In the parameter regime 70 > A/2 the map defined 
by (611 is no longer invertible for all times due to the 


d 1 

Ps{t) = li{i) - Ps{t)], (69) 


existence of zeros of the decoherence function. Thus, the 
model provides an example for a family of maps lying in 
the lower (red) region depicted in Fig. Physically, the 
two-level system reaches the ground state before memory 
effects revive the coherences and the population of the 
upper state. Indeed, also in this case non-Markovianity 
can be traced back to the revival of populations and 
coherences, corresponding to a backflow of information 
from the environment to the system. Since the inverse 
of <!>( does not exist for all times, the criterion based on 
CP-divisibility cannot strictly speaking be applied. This 
is reflected by the fact that the corresponding measure, 


dt 


where Oi with i = x,y,z denote the three Pauli opera¬ 
tors. Considering, e.g., the decoherence dynamics of a 
spin-i particle in a complex environment the rates 7 i(t) 
would correspond to generally time dependent transver¬ 
sal and longitudinal decoherence rates. The dynamical 


map corresponding to (69) can be exactly worked out and 


is given by the random unitary dynamics 


^’t(ps) = XI 


^iPS^i : 


(70) 


i =0 


taking in this case again the expression (591, jumps from 
zero to infinity. 

In the non-resonant case, i.e., for nonzero detuning, 
one can observe a transition from Markovian to non- 
Markovian dynamics for a fixed value of the coupling 
in the weak coupling regime by increasing the detun¬ 
ing A. In this case the map is invertible for all times, 
since the decay rate 7 (f) is finite for all times, so that 
also the measure (48 1 remains well defined. Again, the 


regions of non-Markovianity coincide for all approaches 
since the constraint in order to have CP-divisibility or 
P-divisibility is the same. The pair of states maximiz- 


where (Tq denotes the identity operator and the coeffi¬ 
cients Pi (t) summing up to one are determined from the 
decoherence rates. Introducing the quantities Aij{t) = 
exp[—(ri(t) -l-rj (t))], where as usual we have denoted by 
rfc(t) the time integral of the fc-th decay rate, the coeffi¬ 
cients pi(t) take the explicit formpo,i = (l/ 4 )[l±Ai 2 (t)± 
Ai 3 {t) + A23{t)] and p 2,3 = (l/4)[i T Ai 2 (t) ± Ai3{t) - 
^ 23 (i)]- These equations show that the dynamical map 
actually depends on the sum of the integrals of the de¬ 
cay rates. According to its explicit expression the map 
is completely positive provided the coefficients pi{t) 
remain positive, in which case they can be interpreted as 
a probability distribution and thus characterize the ran¬ 


dom unitary dynamics (70). Note that this can be the 


ing the expression of the measure (341 are however now case even if a decoherence rate stays negative at all times 


given by the projections onto the excited and the ground as in 

Hall et al. ( 

20141 and 

Vacchini et al. 

(2011 

)• 

state. This shows how the optimal pair, which is always As explained in Sec. |IL 2 ^ 

L2 the map is CP-divisible 


given by orthogonal states, does depend on the actual when all of the decoherence rates remain positive for all 
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t ^ 0, so that as soon as at least one of the rates becomes 


negative the measure (48) becomes positive indicating a 


non-Markovian behavior. However, the condition (14) 


for having P-divisibility is weaker. Indeed, in order to 
have P-divisible dynamics only the sum of all pairs of 
distinct decoherence rates has to remain positive, i.e., 
7i(t) + Ijit) ^ 0 for all j ^ i. The same constraint war¬ 
rants monotonicity in time of the behavior of the trace 


distance, so that also in this case the measure (34) of non- 


Markovianity and its generalized version (39) based on 


the Helstrom matrix and corresponding to P-divisibility 
lead to the same result. In the characterization of non- 
Markovianity based on the backflow of information from 
the environment to the system memory effects therefore 
only appear whenever the sum of at least a pair of deco¬ 
herence rates becomes negative. To better understand, 
how this fact can be related to the dynamics of the sys¬ 


tem, one can notice (Chruscihski and Maniscalco 2014) 


that the Bloch vector components {ai{t)), according to 
Eq. (69), obey the equation 




(71) 


where the relaxation times are given by Ti(t) = [jjit) + 
(with all three indices taken to be distinct) and 
correspond to experimentally measurable quantities. Ap¬ 
pearance of non-Markovianity and therefore failure of P- 
divisibility is thus connected to negative relaxation rates, 
corresponding to a rebuild of quantum coherences. For 
the present model one can also easily express the con¬ 
dition leading to a growth of the volume of accessible 
states within the Bloch sphere. In fact, this volume of 
accessible states is proportional to exp[— ri(t)]. As 
a consequence, the dynamics according to the geometric 


criterion by Lorenzo et al. (2013) is Markovian if and only 
if Ell 7i(t) ^ 0 for t ^ 0, which is a strictly weaker cri¬ 
terion with respect to either P or CP-divisibility. 

While in the examples considered above the non- 
Markovianity measure ([M| ) based on the trace distance 
and its generalization ( |39[ ) based on the Helstrom ma¬ 
trix agree in the indication of non-Markovian dynamics, 
a situation in which these criteria are actually different 


has been considered by Chruscihski et al. (2011). 


4. Spin-boson model 


Finally, we briefly discuss a further paradigmatic 
model for dissipative quantum dynamics with many ap¬ 


plications, namely the spin-boson model (Leggett et al. 


1987). The system and the environmental Hamiltonian 


are again given by Eq. (50), while the interaction Hamil¬ 
tonian has the non-rotating-wave structure 


Hi = 




* t 
' 9k^k 


)■ 


(72) 


A^($) 



Q/a;o 


FIG. 5 (Color online) The non-Markovianity measure A/'(<E>) 
defined by Eq. ( |34[ ) for the spin-boson model as a function of 
temperature T and cutoff frequency SI in units of the transi¬ 
tion frequency wo (Clos and Brener 2012). 


Figure]^ shows the non-Markovianity measure M{^) for 
this model obtained from numerical simulations of the 
corresponding second-order time-convolutionless master 
equation, using an Ohmic spectral density of Lorentz- 
Drude shape with cutoff frequency H, reservoir tem¬ 
perature T and a fixed small coupling strength of size 
7 = O.lwo (Clos and Breuer 2012). As can be seen from 
the figure the dynamics is strongly non-Markovian both 
for small cutoff frequencies and for small temperatures 
(note the logarithmic scale of the color bar). The emer¬ 
gence of a Markovian region within the non-Markovian 
regime for small cutoffs and temperatures can be under¬ 
stood in terms of a resonance between the transition fre¬ 
quency and the maximum of an effective, temperature- 
dependent spectral density of the environmental modes. 


B. Applications of quantum non-Markovianity 


In Sec. |HLA| we have illustrated quantum non- 
Markovian behavior by means of simple model systems. 
The goal of the present section is to demonstrate that 
memory effects and their control also open new perspec¬ 
tives for applications. In fact, understanding various as¬ 
pects of non-Markovianity makes it possible to develop 
more refined tools for reservoir engineering and to use 
memory effects as an indicator for the presence of typi¬ 
cal quantum features. Moreover, it becomes possible to 
develop schemes where a small open system is used as 
a quantum probe detecting characteristic properties of 
the complex environment it is interacting with. Here, we 
present a few examples from the recent literature illus¬ 
trating these points. 
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1. Analysis and control of non-Markovian dynamics 

We first consider the model of a spin-^ particle with 
spin operator Sq coupled to a chain of N spin-^ parti¬ 
cles with spin operators s„ (n = 1, 2,..., N), which has 
been investigated by Apollaro et al. (2011). The system 


Hamiltonian is given by Hg = —2hoS§, where Hq denotes 
a local field acting on the system spin Sq. The environ¬ 
mental Hamiltonian represents an AA-Heisenberg spin 
chain with nearest neighbor interactions of strength J in 
a transverse magnetic field h: 


N-l 


N 


He = -2 J («+i + -2h^ 


n—1 


n—1 


The interaction Hamiltonian takes the form 


Hi = -2 Jo (sgs? + . 


(73) 


(74) 


describing an energy-exchange interaction between the 
system spin Sq and the first spin Si of the chain with 
coupling strength Jq. This model leads to a rich dynami¬ 
cal behavior of the system spin and allows to demonstrate 
not only the impact of the system-environment interac¬ 
tion, but also how the open system dynamics changes 
when manipulating the interactions between the con¬ 
stituents of the environment and their local properties. 

It can be shown analytically that optimal state pairs 
for the trace distance measure correspond to antipodal 
points on the equator of the surface of the Bloch sphere. 
Moreover, when the local fields of the system and envi¬ 
ronmental sites are equal, h = ho, one can derive an an¬ 
alytic expression for the rate of change a{t) of the trace 


distance [see Eq. (35l] for optimal state pairs: 


r(t) = -(2A)sgn[J'i(2t)] J 2 ( 2 t), 


(75) 


where sgn is the sign function, and Ji and J ^2 denote 
the Bessel functions of order 1 and 2, respectively. In the 
physically more interesting case h ^ ho the trace distance 


measure [see Eq. (34)] has to be determined numer¬ 

ically. Let us first consider the case, where the system 
local field is zero, ho = 0, the couplings between the spins 
are equal, J = Jg, and we tune the ratio h/J of the local 
field with respect to the spin-spin interaction strength. It 
turns out that the open system is then Markovian only 
at the point h/ J = 1/2, i.e., when the strength of the 
local held is half of the coupling between the spins. Eor 
all other values oi h/J the open system dynamics dis¬ 
plays memory effects with A/’($) > 0. Thus, the point of 
Markovianity separates two regions where memory effects 
are present. This behavior persists for ho ^ 0: Denoting 
the detuning between the local helds hy S = h — ho and 
still keeping J = Jg, the point of Markovianity occurs at 
6h/J = 1/2. On both sides of this point the system again 
exhibits memory effects, but interestingly their proper¬ 
ties are different. Eor Sh/J < 1/2 the system approaches 


a steady state which is independent of the initial state. 
On the other hand, for Sh/J >1/2 the trace distance 
does not decay to zero asymptotically in time. This fea¬ 
ture can be interpreted as information trapping since it 
implies that the distinguishability of states does not van¬ 
ish asymptotically. We also remark that the origin of 
memory effects and the differences in the dynamical be¬ 
havior can be understood in more detail by studying the 


spectrum of the total Hamiltonian (Apollaro et al.. 20111. 


As our second example we discuss applications to a 
physical system which has been studied extensively dur¬ 
ing the last twenty years, namely Bose-Einstein conden¬ 
sates (BECs) which can be manipulated with high preci¬ 
sion by current technologies and thus provide means for 


environment engineering 

Dalfovo et al. 

999 Pethick 

and Smith 2008 fPitaevskii and Stringari 

2003 

). Fol- 

lowing 

Haikka et al. 

(2011 

we consider an impurity-BEC 


system, where the Markovian to non-Markovian transi¬ 
tion can be controlled by tuning the properties of the 
BEC. An impurity atom is trapped within a double-well 
potential, where the left (|A)) and right (|i?)) states rep¬ 
resent the two qubit states. The BEC, which constitutes 
the environment of the qubit, is trapped in a harmonic 
potential. The Hamiltonian for the total impurity-BEC 
system reads 


H = X!^k4ck-b^(^k4 + ^kCk)+X!^^(5k4-f5kCk)- 

k k k 

(76) 

Here, Ak denotes the energy of the Bogoliubov mode Ck 
of the condensate, = |i?)(i?| — \L){L\, while 5k and ^k 
describe the impurity-BEC and the intra-BEC couplings, 
respectively. 

For a background BEC at zero temperature the open 
qubit dynamics is described by a dephasing master equa¬ 


tion of the form of Eq. (54) and it turns out that the 


decoherence rate y{t) can be tuned by changing the 
inter-well distance, and by controlling the dimensional¬ 
ity and the interaction strength (scattering length) of 
the BEC. The detailed studies demonstrate that for a 
3D gas, an increase of the well-separation and of the 
intra-environment interaction leads to an increase of the 
non-Markovianity measure for the impurity dynamics. 
Moreover, the dimensionality of the environment influ¬ 
ences the emergence of memory effects. Studying the 
Markovian to non-Markovian transition as a function of 
the intra-environment interaction strength shows that a 
3D environment is the most sensitive. The Markovian 
- non-Markovian crossover occurs in 3D for weaker val¬ 
ues of the scattering length than in ID and 2D environ¬ 
ments, and a ID environment requires the largest scat¬ 
tering length for the crossover. This behavior can be 
ultimately traced back to the question of how the di¬ 
mensionality of the environment influences the spectral 
density which governs the open system dynamics. In gen¬ 
eral, for small frequencies the spectral density shows a 
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power law behavior, J{ui) cx w®. The spectral density is 
called sub-Ohmic, Ohmic, and super-Ohmic for s < 1, 
s = 1, and s > 1, respectively. For a 3D gas, even with¬ 
out interactions within the gas, the spectral density has 
super-Ohmic character. Introducing interactions within 
the environment then makes this character stronger and 
thereby the 3D gas is most sensitive to the Markovian 
to non-Markovian transition. In contrast, for the ID gas 
with increasing scattering length the spectral density first 
changes from sub-Ohmic to Ohmic and finally to super- 
Ohmic which then allows the appearance of memory ef¬ 
fects. 


2. Open systems as non-Markovian quantum probes 

In addition to detecting, quantifying and controlling 
memory effects, it turns out to be fruitful to ask whether 
the presence or absence of such effects in the dynamics 
of an open quantum system allows to obtain important 
information about characteristic features of a complex 
environment the open system is interacting with. To il¬ 
lustrate this point we present a physical scenario, where 
Markovian behavior indicates the presence of a quan¬ 
tum phase transition in a spin environment ( Sachde'^ 
2011). Consider a system qubit, with states |g) and |e), 
which interacts with an environment described by a one¬ 
dimensional Ising model in a transverse field (iQuan et al.\ 


2006). The Hamiltonians of the environment and the 


qubit-environment interaction are given by 


= ( 77 ) 

3 

Hj = -J5\e){e\Y,^l (78) 


where J is a microscopic energy scale, A a dimensionless 
coupling constant describing the strength of the trans¬ 
verse field, and are Pauli spin operators. The sys¬ 
tem qubit is coupled with strength 5 to all environmental 
spins. The model yields pure dephasing dynamics de¬ 
scribed by a decoherence function Git) which, for a pure 
environmental initial state |<I>), is given by 


iHnt ^ — 


Git) = 




(79) 


where the effective environmental Hamiltonians are 
Ha = He + (a| Hj |a) with a = g,e. We note that ac¬ 


cording to Eq. (79) the decoherence function is directly 
linked to the concept of the Loschmidt echo L{t) (Cue- 
chietti et al. 2003) which characterizes how the environ 


ment responds to perturbations by the system ( 

Gucchi- 

etti et aLl|2004[ Gorin et al.l|20061|Goussev et al.\ 2008 

Peres 

1984 

). In fact, we have L{t) = G(t)p* and, hence, 


there is also a direct connection between the Loschmidt 
echo and the evolution of the trace distance for optimal 
state pairs, i.e. D{pgit), pg{t)) = ^JLit). 



FIG. 6 (Co lor online) The non-Markovianity measure A/'(‘I>) 
of Eq. (341 for a qubit coupled to a one-dimensional Ising 
chain in a transverse field as a function of field strength A* 
and particle number N. The measure vanishes along the black 
line and is nonzero everywhere else (Haikka et al. 20121. 


The system-environment interaction (78) leads to the 
effective transverse field A* = A -I- d when the system is 
in the state |e). At the critical point of the spin environ¬ 
ment (A* = 1) the spin-spin interaction and the external 
field are of the same size and the environment is most 
sensitive to external perturbations. This is reflected in a 
quite remarkable manner in the dynamics of the system 
qubit as can be seen from Fig. ( [Haikka et al. 2012). 
At the point A* = 1 the system qubit experiences strong 
decoherence. As a matter of fact, it turns out that the dy¬ 
namics of the system qubit always shows non-Markovian 
behavior except at the critical point of the environment 
where the dynamics is Markovian. Therefore, the sys¬ 
tem qubit acts as a probe for the spin environment and 
Markovian dynamics represents a reliable indicator of en¬ 
vironment criticality. It is important to note that even 
though, strictly speaking, quantum phase transitions re¬ 
quire to take the thermodynamic limit, the above results 
hold for any number of environmental spins. 


IV. IMPACT OF CORRELATIONS AND EXPERIMENTAL 
REALIZATIONS 


In this section we study the effect of two fundamental 
kinds of correlations on the non-Markovian behavior of 
open systems, namely correlations between the open sys¬ 
tem and its environment and correlations within a com¬ 
posite environment. We conclude the section by giving 
an overview of the experimental status of non-Markovian 
open systems. 
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A. Initial system-environment correlations 

1. Initial correlations and dynamical maps 

Up to this point we have discussed the theory of open 
quantum systems via the concept of a dynamical map, 
which in Sec. III. ATT] we wrote down in terms of the initial 
state of the environment and of the total Hamiltonian. 
In doing so we assumed that the open system and the 
environment are independent at the initial time of prepa¬ 
ration [see Eqs. Q and 0]. However, experimentalists 
often study fast processes in strongly coupled systems, 
where this independence is rarely the case, as originally 


pointed out by Pechukas (1994). Going beyond this ap¬ 


proximation brings us up against the problem of describ¬ 
ing the dynamical properties of open quantum systems 
without dynamical maps. 

The theory of completely positive maps plays a crucial 
role in many fields of quantum physics and thus a gen¬ 
eral map description for open systems, even in the pres¬ 
ence of initial correlations, would be very useful. Thus, 
naturally, a lot of effort has been put in trying to ex¬ 


tend the formalism along these lines (Alicki 1995 Sha- 


bani and Lidar 

2009 

Stelmachovic and Buzek 

2001 


In this framework the aim has been to determine un¬ 
der which conditions the system dynamics can be de¬ 
scribed by means of completely positive maps if initial 
correlations are present. Unfortunately, no generally ac¬ 
knowledged answer to this question has been found and 


the topic is still under discussion ( 

Brodutch et al. 

2013 

Rodriguez-Rosario et al. 

2010 

Shabani and Lidar 

o 

o 


As useful as a description in terms of completely pos¬ 
itive maps would be, one may as well ask the question: 
How do the initial correlations influence the dynamics? 
Especially, since an experimentalist is often restricted to 
looking only at a subset of system states, he would want 
to conclude something about the system based on the 
dynamics of this subset. It is not possible to construct a 
map, but perhaps something else could be concluded. In 
the following we present an approach, in which the dy¬ 
namics of information flow between the system and the 
environment is studied. We will see that the initial cor¬ 
relations modify the information flow and consequently 
can be witnessed from the dynamical features of the open 
system. 

Let us have a closer look at the dynamics of the infor¬ 


mation inside the open system lint (t) , defined in Eq. (301, 


for a general open system, where initial correlations could 
be present. Suppose there are two possible preparations 
giving rise to distinct system-environment initial states 
Pse(^) Since the total system and the envi¬ 

ronment evolve under unitary dynamics, we have 

— 2iint(0) = lext(O) — lext(i) ^ 2iext(0), (80) 

where Text {t) is the information outside the open system 


defined in Eq. (31). This inequality, as simple as it is. 


reveals one very important point: The information in the 
open system lint(i) can increase above its initial value 
2^int(0) only if there is some information initially outside 
the system, i.e. lext(O) > 0. Obviously, the contraction 
property (28) for completely positive maps is a special 
case of the inequality (80) which occurs when the system 
and the environment are initially uncorrelated and the 
environmental state is not influenced by the system state 


preparation, i.e. = pg ®pE- This leads toXcxt(O) = 
0 and, thus, Eq. (fSOj) reduces to Eq. (p9). 


2. Local detection of initial correlations 


We found that initial information outside the open sys¬ 
tem can lead to an increase of trace distance. But how 
can this be used to develop experimental methods to de¬ 
tect correlations in some unknown initial state P 5 £? To 
this end, let us combine (801 and the inequality (33) for 
the initial time in order to reveal the role of initial cor¬ 
relations more explicitly (Laine et al.\ 2010b): 


iZint(0 — 2^int(0) ^ i4(p5^(0), p;g(0) ® Pb(0)) 

+-O(Ps'£;(0)) Ps(0) G p|;(0)) 
+Z?(pi;(0),p|(0)). (81) 


This inequality clearly shows that an increase of the trace 
distance of the reduced states implies that there are ini¬ 
tial correlations in P 5 £;( 0 ) or p|£;(0), or that the initial 
environmental states are different. 



FIG. 7 (Color online) Schematic picture of the behavior of 
the trace distance with and without initial correlations. If 
the trace distance between two states, prepared such that 
pIb(O) = (^®7)psb( 0), never exceeds its initial value (dashed 
black line), the dynamics does not witness initial correlations. 
If, on the other hand, the trace distance for such pair of states 
increases above the initial value (solid red line), we can con¬ 
clude the presence of initial correlations. 


Let us now assume that one can perform a state to¬ 
mography on the open system at the initial time zero 
and at some later time t, to determine the reduced states 
Pg(0) and Pg{t). In order to apply inequality (81) to 
detect initial correlations we need a second reference 
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state pKO), which has the same environmental state, i.e. 
p%{0) = Pe{0)- This can be achieved by performing a 
focal quantum operation on pg^(O) to obtain the state 


Ps'_e( 0) — (A (g) /)pg£;(0). 


(82) 


The operation A acts locally on the variables of the open 
system, and may be realized, for instance, by the mea¬ 
surement of an observable of the open system, or by a 
unitary transformation induced, e.g., through an exter¬ 
nal control field. Now, if Pg^i^) i® uncorrelated then also 
p|^( 0 ) is uncorrelated since it has been obtained through 
a local operation. Therefore, for an initially uncorre¬ 
lated state the trace distance cannot increase according 
to inequality 0 . Thus, any such increase represents a 
witness for correlations in the initial state P 5 £:( 0 ), as is 
illustrated in Fig. We note that this method for the 
focal detection of initial correlations requires only focal 
control and measurements of the open quantum system, 
which makes it feasible experimentally and thus attrac¬ 
tive for applications. In fact, experimental realizations of 


the scheme have been reported recently, see Sec. IV.C.2 


Furthermore, Smirne et al. (2010) have studied how the 


scheme can be employed for detecting correlations in 
thermal equilibrium states. 




1 tTEoU{t) ^ 

-Pse(O) --—--► Ps(t) 


A(gil 


Pl( 0 ) ◄- 


trs 


Pse(O) ■ 


tlEo!7(t) ^ 

—--► Ps(^) 


FIG. 8 (Color online) Scheme for the local detection of 
system-environment quantum correlations in some state 
Ps£;(0)i employing a local dephasing operation A to gener¬ 
ate a second reference state PeeW- 

The above strategy can even be used in order to lo¬ 
cally detect a specific type of quantum correlations, i.e. 
system-environment states with nonzero quantum dis¬ 


cord (Modi et al. 2012). This is achieved by taking the 


focal operation A in Eq. (82) to be a dephasing operation 


leading to c omplete decoherence in the eigenb asis of the 
state ^ 5 ( 0 ) ( [Gessner and Brener 2011, 2013). The ap¬ 
plication of this dephasing operation destroys all quan¬ 
tum correlations of p; 5 ^( 0 ), while leaving invariant its 
marginal states ^ 5 ( 0 ) and p\{G) (see Fig. |^. Thus, one 
can use the quantity C(p 5 ^( 0 )) = D{pg^[Q), pgE:{^)) 9 , 

measure for the quantum correlations in the initial state 
Pse{^) (Luo 2008). Since the trace distance is invariant 
under unitary transformations and since the partial trace 
is a positive map, the corresponding local trace distance 
provides a lower bound of C(/ 9 ;g^( 0 )) for 
all times t ^ 0. Hence, we have 


maxD{pg{t),psit)) ^ CipsEiO))- 


(83) 


The left-hand side of this inequality yields a locally acces¬ 
sible lower bound for a measure of the quantum discord. 
Experimental realizations of the scheme are briefly de¬ 
scribed in Sec. EY Most recent studies suggest that 
the method could also be applied for detecting a criti¬ 
cal point of a quantum phase transition via dynamical 


monitoring of a single spin alone (Gessner et al. 2014a). 


B. Nonlocal memory effects 

We now return to study the case of initially uncor¬ 
related system-environment states, for which the usual 
description in terms of dynamical maps can be used. So 
far, we have concentrated on examples with a single sys¬ 
tem embedded in an environment. However, many rel¬ 
evant examples in quantum information processing in¬ 
clude multipartite systems for which nonlocal proper¬ 
ties such as entanglement become important and it is 
therefore essential to study how non-Markovianity is in¬ 
fluenced by scaling up the number of particles. 

The question of additivity of memory effects with re¬ 
spect to particle number was recently raised by |Addis| 
et al. (2013, 2014); and Fanchini et al. (2013) and it was 


found that the different measures for non-Markovianity 
have very distinct additivity properties. Indeed, the re¬ 
search on multipartite open quantum systems is yet in its 
infancy and a conclusive analysis remains undone. In the 
following, we will discuss a particular feature of bipartite 
open systems: the appearance of global memory effects 
in the absence of local non-Markovian dynamics. This 
is at variance with the standard situation in which the 
enlargement of considered degrees of freedom leads from 
a non-Markovian to a Markovian dynamics (Martinazzo 


et al. 2011 ). 


We will now take a closer look at nonlocal maps, for 
which memory effects may occur even in the absence of 


local non-Markovian effects. In Laine et al. (2012 2013) 


it is shown that such maps may be generated from a 
local interaction, when correlations between the environ¬ 
ments are present and that they may exhibit dynamics 
with strong global memory effects although the local dy¬ 
namics is Markovian. The nonlocal memory effects are 
studied in two dephasing models: in a generic model of 


qubits interacting with correlated multimode fields ( 

Wifi- 

mann and Breuer 

2013 

Wifimann and Breuer 

2014 

) and 


in an experimentally realizable model of down converted 
photons traveling through quartz plates, which we will 
discuss later in detail. 

Before we turn to discuss any experimental endeavor to 
detect non-Markovian dynamics, let us discuss in more 
rigor about the generation and possible applications of 
nonlocal memory effects. Consider a generic scenario, 
where there are two systems, labeled with indices i = 1 , 2 , 
which interact locally with their respective environments. 
The dynamics of the two systems can be described via 
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the dynamical map 


C. Experiments on non-Markovianity and correlations 


pf{t) = <^>rpf{0) 


= tiE 




where U^^{t) = Ug^{t) ®Ug^{t) with Ug^{t) describing 
the local interaction between the system i and its environ¬ 
ment. If the initial environment state p^{0) factorizes, 
i.e. Pb (0) = Pe{0) ®/o|;(0), also the map factorizes. 
Thus, the dynamics of the two systems is given by a lo¬ 
cal map ® $(. On the other hand, if pjj (0) 

exhibits correlations, the map cannot, in general, be 
factorized. Consequently, the environmental correlations 
may give rise to a nonlocal process even though the in¬ 
teraction Hamiltonian is purely local (see Fig. |^. 


Up to this point we have not yet discussed any ex¬ 
perimental aspects of the detection of memory effects. 
In the framework of open quantum systems the environ¬ 
ment is in general composed of many degrees of freedom 
and is therefore difficult to access or control. To perform 
experiments on systems where the environment induced 
dynamical features can be controlled is thus challenging. 
However, in the past years clever schemes for modifying 
the environment have been developed allowing the estab¬ 
lishment of robust designs for noise engineering. In this 
section we will briefly review some of the experimental 
platforms where a high level of control over the envi¬ 
ronment degrees of freedom has been accomplished and 
non-Markovian dynamics observed and quantified. 



FIG. 9 (Color online) Schematic picture of a system with 
nonlocal memory effects. The open systems 1 and 2 locally 
interact with their respective environments. Initial correla¬ 
tions between the local environments cause the occurrence of 
nonlocal memory effects. 

For a local dynamical process, all the dynamical prop¬ 
erties of the subsystems are inherited by the global sys¬ 
tem, but naturally for a nonlocal process the global dy¬ 
namics can display characteristics absent in the dynamics 
of the local constituents. Especially, for a nonlocal pro¬ 
cess, even if the subsystems undergo a Markovian evo¬ 
lution, the global dynamics can nevertheless be highly 
non-Markovian as was shown for the dephasing models 


by Laine et al. (2012 2013). 


Thus, a system can globally recover its earlier lost 
quantum properties although the constituent parts are 
undergoing decoherence and this way initial environmen¬ 
tal correlations can diminish the otherwise destructive 
effects of decoherence. This feature has been further 
deployed in noisy quantum information protocols, such 


1. Control and quantification of memory effects in photonic 
systems 

Quantum optical experiments have for many decades 
been the bedrock for testing fundamental paradigms of 
quantum mechanics. The appeal for using photonic sys¬ 
tems arises from the extremely high level of control allow¬ 
ing, for example, controlled interactions between different 
degrees of freedom, preparation of arbitrary polarization 
states and a full state tomography. Needless to say, pho¬ 
tons thus offer an attractive experimental platform also 
for studying non-Markovian effects. 


Liu et al. (2011) introduced an all-optical experiment 


which allows through careful manipulation of the initial 
environmental states to drive the open system dynam¬ 
ics from the Markovian to the non-Markovian regime, 
to control the information flow between the system and 
the environment, and to determine the degree of non- 
Markovianity. In the experiment the photon polariza¬ 
tion degree of freedom (with basis states \H) and \ V) for 
horizontal and vertical polarization, respectively) plays 
the role of the open system. The environment is repre¬ 
sented by the frequency degree of freedom of the pho¬ 
ton (with basis {|w)}^^g), which is coupled to the sys¬ 
tem via an interaction induced by a birefringent mate¬ 
rial (quartz plate). The interaction between the polar¬ 
ization and frequency degrees of freedom in the quartz 
plate of thickness L is described by the unitary operator 
U{t) I A) 0 Iw) = I A) 0 |w), where n\ is the refrac¬ 

tion index for a photon with polarization X = H,V and 
t = L/c is the interaction time with the speed of light c. 
The photon is initially prepared in the state \ip) 0 |x), 
with Itp) = a \H) + P \V) and |x) = f dujf{oj) |a;), where 
/(w) gives the amplitude for the photon to be in a 


as teleportation ( 

Laine et al. 

2014 

) and entanglement mode with frequency uj. The quartz plate leads to a 

distribution (IXiang et al. 

201 

4), suggesting that non- pure decoherence dynamics of superpositions of polariza- 


Markovianity could be a resource for quantum informa¬ 
tion. 


decoherence function G{t) = J dw\f{uj)pe 


2^iu:Ant 
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An = ny — uh- An optimal pair of states maximizing 
the non-Markovianity measure (341 for this map is given 
by |V^i, 2 ) = ± |14)) and the corresponding trace 

distance is D{pg{t), p‘g{t)) = |G'(t)|. Clearly, controlling 
the frequency spectrum |/(a;)p changes the dynamical 
features of the map. 




FIG. 10 (Color online) Non-Markovian dynamics arising from 
an engineered frequency spectrum, (a) The frequency spec¬ 
trum of the initial state for various values of the tilting angle 6 
of the cavity, (b) The change of the trace distance as functions 
of the tilting angle 6. The transition from the non-Markovian 
to the Markovian regime occurs at 9 « 4.1°, and from the 
Markovian to the non-Markovian regime at 6 « 8.0°, corre¬ 
sponding to the occurrence of a double peak structure in the 
frequency spectrum (a) (Liu et al. 20111. 


The modification of the photon frequency spectrum is 
realized by means of a Fabry-Perot cavity mounted on a 
rotator which can be tilted in the horizontal plane. The 
shape of the frequency spectrum is changed by chang¬ 
ing the tilting angle 9 of the Fabry-Perot cavity, as can 
be seen in Fig. [TO^ . Further, the frequency spectrum 
changes the dephasing dynamics such that the open sys¬ 
tem exhibits a reversed flow of information and thus 
allows to tune the dynamics from Markovian to non- 
Markovian regime (see Fig. [ToIj). In the experiment, full 


state tomography can further be performed thus allowing 
a rigorous quantification of the memory effects. 

Also, a series of other experimental studies on non- 
Markovian dynamics in photonic systems have been per¬ 
formed recently. Tang et al. (2012) report a measure¬ 
ment of the non-Markovianity of a process with tunable 
system- environment interaction, Cialdi et al. (2011) ob¬ 
serve controllable entanglement oscillations in an effec¬ 


tive non-Markovian channel and in (Chiuri et al.. 2012 


Jin et al. 20141 simulation platforms for a wide class of 


non-Markovian channels are presented. 


2. Experiments on the local detection of correlations 


A number of experiments have been carried out in or¬ 
der to demonstrate the local scheme for the detection 
of correlations between an open system and its environ¬ 
ment described in Sec. IIV.A.21 Photonic realizations of 


this scheme have been reported in (Li et aZ. 2011 Smirne 


et al. 2011), where the presence of initial correlations be¬ 


tween the polarization and the spatial degrees of freedom 
of photons is shown by the observation of an increase of 
the local trace distance between a pair of initial state. 

As explained in Sec. |IV.AT2l it is also possible to reveal 
locally the presence of quantum correlations represented 
by system-environment states with nonzero discord if the 


local operation A in Eq. (821 is taken to induce complete 


decoherence in the eigenbasis of the open system state. 



the discrimination between quantum and classical corre¬ 
lations, and applied this extension to a photonic experi¬ 
mental realization. 

All experiments mentioned so far employ photonic de¬ 
grees of freedom to demonstrate non-Markovianity and 
system-environment correlations. The first experiment 
showing these phenomena for matter degrees of freedom 


has been described by Gessner et al. (2014b). In this ex¬ 


periment nonclassical correlations between the internal 
electronic degrees of freedom and the external motional 
degrees of freedom of a trapped ion have been observed 


and quantified by use of Eq. (83). Important features of 


the experiment are that the lower bounds obtained from 
the experimental data are remarkably close to the true 
quantum correlations present in the initial state, and that 
it also allows the study of the temperature dependence 
of the effect. 


3. Non-Markovian quantum probes detecting nonlocal 
correlations in composite environments 

In Sec. lIV.Bl we demonstrated that initial correlations 
between local parts of the environment can lead to non- 
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local memory effects. In Liu et al. (2013a I such nonlocal 


memory effects were experimentally realized in a pho¬ 
tonics system, where manipulating the correlations of the 
photonic environments led to non-Markovian dynamics of 
the open system. The experimental scheme further pro¬ 
vided a controllable diagnostic tool for the quantification 
of these correlations by repeated tomographic measure¬ 
ments of the polarization. 

Let us take a closer look at the system under study in 
the experiment. A general pure initial polarization state 
of a photon pair can be written as \ 1 p 12 ) = a\HH) + 
b\HV) + c\VH) + d\VV) and all initial states of the 
polarization plus the frequency degrees of freedom are 
product states 

1^(0)) = I'012) ® J duJiduj2 9{i^i,i^2)\i^i,i^2) , (85) 


where g{uJi,u} 2 ) is the probability amplitude for photon 1 
to have frequency Wi and for photon 2 to have frequency 
W 2 , with the corresponding joint probability distribution 
P{oji,uj 2 ) = | 5 (wi,W 2 )p. If the photons pass through 
quartz plates, the dynamics can be described by a gen¬ 
eral two-qubit dephasing map, where the different de¬ 
coherence functions can be expressed in terms of Fourier 
transforms of the joint probability distribution P(a;i, a; 2 ). 
For a Gaussian joint frequency distribution with identical 
single frequency variances C and correlation coefficient 
K, the time evolution of the trace distance correspond¬ 
ing to the Bell-state pair = -^{\HH) ± \VV)) is 
found to be 


Dit) = exp 




( 86 ) 


where ti denotes the time photon i has interacted with 
its quartz plate up to the actual observation time t. For 
uncorrelated photon frequencies we have K = 0 and the 
trace distance decreases monotonically, corresponding to 
Markovian dynamics. However, as soon as the frequen¬ 
cies are anticorrelated, K < 0, the trace distance is non¬ 
monotonic which signifies quantum memory effects and 
non-Markovian behavior. On the other hand, the local 
frequency distributions are Gaussian and thus for the sin¬ 
gle photons the trace distance monotonically decreases. 
Therefore we can conclude, that the system is locally 
Markovian but globally displays nonlocal memory effects. 

In the experiment two quartz plates act consecutively 
for the photons, and the magnitude of the initial anti¬ 
correlations between the local reservoirs is tuned. After 
the photon exits the quartz plates, full two-photon polar¬ 
ization state tomography is performed and by changing 
the quartz plate thicknesses, the trace distance dynam¬ 
ics is recovered. From the dynamics it is evident that 
initial environmental correlations influence the quantum 
non- M arkovianity. 

A further important aspect of the experimental scheme 
is that it enables to determine the frequency correlation 


coefficient K of the photon pairs from measurements per¬ 
formed on the polarization degree of freedom. Thus by 
performing tomography on a small system we can ob¬ 
tain information on frequency correlations, difficult to 
measure directly. Thus, in this context, the open system 
(polarization degrees of freedom) can serve as a quantum 
probe which allows us to gain nontrivial information on 
the correlations in the environment (frequency degrees of 
freedom). 


V. SUMMARY AND OUTLOOK 

In this Golloquium we have presented recent advances 
in the definition and characterization of non-Markovian 
dynamics for an open quantum system. While in the clas¬ 
sical case the very definition of Markovian stochastic pro¬ 
cess can be explicitly given in terms of constraints on the 
conditional probabilities of the process, in the quantum 
realm the peculiar role of measurements prevents a direct 
formulation along the same path and new approaches are 
called for in order to describe memory effects. We have 
therefore introduced a notion of non-Markovianity of a 
quantum dynamical map based on the behavior in time of 
the distinguishability of different system states, as quan¬ 
tified by their trace distance, which provides an intrinsic 
characterization of the dynamics. This approach allows 
to connect non-Markovianity with the flow of informa¬ 
tion from the environment back to the system and nat¬ 
urally leads to the introduction of quantum information 
concepts. It is further shown to be connected to other 
approaches recently presented in the literature building 
on divisibility properties of the quantum dynamical map. 
In particular, a generalization of the trace distance cri¬ 
terion allows to identify Markovian time evolutions with 
quantum evolutions which are P-divisible, thus leading 
to a clear-cut connection between memory effects in the 
classical and quantum regimes. As a crucial feature the 
trace distance approach to non-Markovianity can be ex¬ 
perimentally tested and allows for the study of system- 
environment correlations, as well as nonlocal memory ef¬ 
fects emerging from correlations within the environment. 

We have illustrated the basic feature of the considered 
theoretical approaches to quantum non-Markovianity, 
typically leading to a recovery of quantum coherence 
properties, by analyzing in detail simple paradigmatic 
model systems, and further discussing more complex 
models which show important connections between non- 
Markovianity of the open system dynamics and features 
of the environment. It is indeed possible to obtain infor¬ 
mation on complex quantum systems via study of the dy¬ 
namics of a small quantum probe. The study of the time 
development of the trace distance between system states 
can be shown to provide an indication of the presence 
of initial correlations between system and environment, 
thus providing a powerful tool for the local detection of 
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correlations. Further considering suitable local opera¬ 
tions, such correlations can be further characterized and, 
in particular, one can distinguish quantum and classical 
correlations. 

Experimentally, it is in general obviously very difficult 
to control in detail the environmental degrees of freedom 
and therefore experiments on non-Markovian quantum 
systems are still in a state of infancy. Still, some so¬ 
phisticated schemes for modifying the environment have 
been developed, thus leading to the realization of first 
proof-of-principle experiments in which to test the study 
of non-Markovianity, its connection with crucial features 
of the environment, as well as the capability to unveil 
system-environment correlations by means of local ob¬ 
servations on the system. 

In recent years non-Markovian quantum systems have 
been enjoying much attention both due to fundamental 
reasons and foreseeable applications, as shown by the 
rapid growth in the related literature. Indeed, the po¬ 
tential relevance of memory effects in the field of com¬ 
plex quantum systems and quantum information has led 
to an intense study, and we have pointed to some high¬ 
lights in this novel, so far fairly unexplored field of non- 
Markovianity. However, there are numerous open ques¬ 
tions yet to be studied. The latter include fundamental 
questions such as the mathematical structure of the space 
of non-Markovian quantum dynamical maps, the role of 
complexity in the emergence of memory effects or the 
relevance of non-Markovianity in the study of the border 
between classical and quantum aspects of nature, as well 
as more applied issues like the identification of the en¬ 
vironmental features or system-environment correlations 
which can indeed be detected by means of local observa¬ 
tions on the system. 

The theoretical and experimental investigations of 
non-Markovian quantum systems outlined in this Col¬ 
loquium pave the way for new lines of research by both 
shedding light on fundamental questions of open quan¬ 
tum systems as well as by suggesting novel applications 
in quantum information and probing of complex systems. 
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